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Abstract Radio synchrotron emission, its polarization and Faraday rotation of the 
polarization angle are powerful tools to study the strength and structure of magnetic 
fields in galaxies. Unpolarized synchrotron emission traces isotropic turbulent fields 
which are strongest in spiral arms and bars (20 — 30 /rG) and in central starburst re¬ 
gions (50— 100/rG). Such fields are dynamically important; they affect gas flows and 
drive gas inflows in central regions. Polarized emission traces ordered fields, which 
can be regular or anisotropic turbulent, where the latter originates from isotropic 
turbulent fields by the action of compression or shear. The strongest ordered fields 
(10—15 /rG) are generally found in interarm regions. In galaxies with strong density 
waves, ordered fields are also observed at the inner edges of spiral arms. Ordered 
fields with spiral patterns exist in grand-design, barred and flocculent galaxies, and 
in central regions. Ordered fields in interacting galaxies have asymmetric distribu¬ 
tions and are a tracer of past interactions between galaxies or with the intergalactic 
medium. - Faraday rotation measures (RM) of the diffuse polarized radio emission 
from galaxy disks reveal large-scale spiral patterns that can be described by the super¬ 
position of azimuthal modes; these are signatures of regular fields generated by mean- 
field dynamos. “Magnetic arms” between gaseous spiral arms may also be products 
of dynamo action, but need a stable spiral pattern to develop. Helically twisted field 
loops winding around spiral arms were found in two galaxies so far. Large-scale field 
reversals, like the one found in the Milky Way, could not yet be detected in external 
galaxies. In radio halos around edge-on galaxies, ordered magnetic fields with X- 
shaped patterns are observed. - The origin and evolution of cosmic magnetic fields, 
in particular their first occurrence in young galaxies and their dynamical importance 
during galaxy evolution, will be studied with forthcoming radio telescopes like the 
Square Kilometre Array. 
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1 The role of magnetic fields in spiral galaxies 


Whether and how magnetic fields influence the formation and evolution of galaxies is 
still unknown. Most numerical models of the interstellar medium in galaxies ignored 
magnetic fields. Magneto-hydrodynamical (MHD) models of evolving galaxies are 


computationally challenging and require various simplifications (|Wang & Abel |2009 
Pakmor & Springel|2013[ Pakmor et al.|2014| l. 


Much more is known about magnetic fields in the interstellar medium (ISM) of 
nearby galaxies and the Milky Way, thanks to modern radio telescopes, the construc¬ 
tion of sensitive and wide-hand receiving systems and the development of new meth¬ 
ods of observation and data analysis in the radio range. Major progress was also 


achieved in numerical modelling of the magnetized ISM (e.g. de Avillez & Bre- 

itschwerdt 

2005 

Gissinger et al.||2009| Hill et al.||2012 Gressel et al.||2013 Gent 

et al.|20I3 

Mac 

lida et al.|20I3|l, using a variety of codes and assumptions. Not sur- 


prisingly, the results differ and are not always consistent with observations, which 
calls for further efforts on both sides. 

A widespread argument to neglect galactic magnetic fields in the past was based 
on assuming field strengths of a few /rG taken from Faraday rotation data in our 
Milky Way, which trace only a small component of the total field. The dynamical im¬ 
portance of galactic fields, with typical strengths of 15 fiG, is now widely appreciated 
in many fields of astrophysics. Magnetic and cosmic-ray pressures control the over¬ 
all star-formation rate ( Birnboim et al.||2 015|l. Magnetic fields are im portant for the 
gas dynamics of molecular clouds ( Planck Collaboration et al.|2015a I. With stronger 
fields fewer cloud cores are formed but with larger masses (Vazquez-Semadeni et al. 
2005 Price & Bate 2008| l. Magnetic fields control the density and propagation of 
cosmic rays. Together with cosmic rays they can provide the pressure to drive fast 
outflows of hot gas, in particular in galaxies with high star-formation rates in the 
early Universe ( |Hanasz et al.||20l3] l. Cosmic rays propagating along the field lines 
in galaxy halos may heat the warm ionized gas and explain the optical line ratios 
( Wiener et al.|2013] l. Outflows from starburst galaxies may have magnetized the in- 
tergalactic medium ( |Bertone et al.|200^ . 

The detection of ultrahigh-energy cosmic rays (UHECRs) reaching the Earth with 
the AUGER observatory and the possibly anisotropic distribution of their arrival di¬ 
rections (Aab et al. 2015| l calls for a proper model of particle propagation. As UHECR 
particles are deflected by large-scale regular fields and scattered by turbulent fields, 
the structure and the extent of the fields in the disk and halo of the Milky Way need to 
be known, but are hard to measure from our inside location. The view onto external 
spiral galaxies can help. 


2 Origin of galactic magnetic fields 

Building up magnetic fields in galaxies is a three-stage process: (1) seeding, (2) am¬ 
plifying, (3) ordering and sustaining. 

(1) Seed fields can be “primordial”, generated in the early Universe ( |Durrer &| 
|Neronov||2013|l, or they may originate from later epochs, e.g. during cosmological 
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structure formation by the Weibel instability ( LazaretalJ2009jfffom plasma fluctua¬ 
tions in protogalaxies (Schlickeiser 2012[|Schlickeiser & Felten 20131, from injection 
by the first stars ( |Bisnovatyi-Kogan et al.||1973 1 or jets generated by the first black 
holes (|Rees|2005|l or from the Biermann mechanism in the first supernova remnants 


( Hanayama et al.|2005| l. 

Indication for a weak all-pervasive intergalactic seed field as a relic from the early 
Universe comes from the non-detection of GeV 7 -ray emission with the FERMI 
satellite from blazars, which were observed at TeV energies with the HESS obser¬ 
vatory. The secondary particles generated by the TeV photons may be defle cted by 
intergalactic fields of at least 10^ G strength and a high volume filling factor (Dolag 
et al.| 2011 | l. However, fluctuations in the intergalactic plasma may also disperse the 
7 -ray emission (Broderick et al.|2012|l (see also the discussion in Durrer & Neronov 
( 2013| l). Analysis of the CMB power spectra data from PLANCK gave an upper limit 
of about 5 nG on a (comoving) scale of 1 Mpc ( Planck Collaboration et al.|2015bl l. 

(2) An efficient source of field amplification is turbulence in the gas driven by 
supernova explosions ( Ferriere|1996| l or by spiral shocks ( Kim et al.|2006| l, called the 
small-scale dynamo. In protogalaxies this mechanism can amplify weak seed fields to 
several pG strength (the energy level of turbulence) within less than 10^ yr (Kulsrud 


et al.|I99'7} Schleicher et al.|20I0 Beck et al.|20I2a Rieder & Teyssier|20i5| l. The 

resulting field is turbulent. 

(3) The final and most time consuming stage is the ordering of the turbulent field. 
Numerical models of evolving galaxies show fast field amplification and ordering 
with help of differential rotation ( |Wang & Abel|2009] [Kotarba et al. 2009|l, possibly 
supported by the magneto-rotational instability (MRI) (Gressel et a LI20131 |Pakmor| 


& Springel|2013 1. In these models the ordered magnetic field forms spiral arm seg¬ 
ments, but it has frequent reversals in azimuthal and radial directions and hence is 
anisotropic turbulent (Sect. HI}- 

The most promising mechanism to sustain magnetic fields and generate large- 
scale regular fields from turbulent fields in the interstellar medium of galaxies is the 
a — Gl dynamo ( [Beck et al.|1996| l. It is based on differential rotation (X2), expanding 
gas flows, driven by supernova explosions ( [Ferriere 1996 Gressel et al.|2013] l or cos¬ 
mic rays ( [Hanasz et al. 20091, carrying magnetic fields that are twisted by the Coriolis 
force (a-efifect), and on magnetic diffusivity (rj) driven by turbulence. The “mean- 
field” approximation of the a — dynamo equation allows analytical solutions by 
mathematically separating the large-scale and small-scale parts of the velocity held 
and the magnetic held. The “mean-field” dynamo equation for the regular held B 


= V X (v X B) -I- V X aB +T] V^B (1) 

describes the held amplification by the induction term V x (v x B), where v is the 
large-scale velocity, the gain term V x aB and the loss term rj V^B. MHD model¬ 
ing with high spatial resolution showed that the scale separation of the mean-field 
approximation is reasonable ( [Gent et al.|2013| l. 

Solutions of the mean-field dynamo equation are described by modes m with dif¬ 
ferent azimuthal symmetries in the disk plane (axisymmetric, bisymmetric, quadri- 
symmetric, etc.) and two different vertical symmetries (even parity or odd parity) 
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perpendicular to the disk plane. Several modes can be excited in the same object, 
where lower modes have shorter growth times and larger amplitudes. In flat, rotat¬ 
ing objects like galaxy disks, the strongest mode is composed of a toroidal held of 
axisymmetric spiral shape within the disk (m = 0), without sign reversals across the 
equatorial plane, and a weaker poloidal held of even symmetry. Azimuthal dynamo 
modes can be identihed observationally from the patterns of polarization angles and 
Faraday rotation measures of the diffuse polarized emission from galaxy disks (see 
Sect. [ 43 . 

Mean-held models of the a — dynamo in galaxy disks predict that turbulent 


helds are transformed into large-scale regular helds within a few 10^ yr (Beck et al. 


1994[ Arshakian et al.|[200^ Rodrigues et al.|20T5] l. Field reversals from the early 
phases may survive until today (|Moss et al.|2012[). Global numerical models of galax¬ 


ies (Gissinger et al.|[2()09[|Hanasz et al.||2009 Kulpa-Dybel et al.|2011[|Siejkowski| 


|et al. 12014 1 conhrmed the basic results of the mean-held models. 

The twisted held loops needed for the operation of the a — dynam o generate 
large-scale helicity.lHAs total helicity in a system is a conserved quantity (Branden- 
burg & Subramanianl 2005| l, small-scale helds with opposite helicity are generated 
which quench dynamo action, unless they are removed from the system (|Shukurovl 


et al.||20()6|. Flence, moderately fast outhows are essential for effective a — Q dy¬ 


namo action, but strong outhows can suppress dynamos. This leaves a certain range 
of outhow velocities for optimal dynamo action ( Rodrigues et al.|2015 1. As outhows 
are strongest above spiral arms, the a — Q dynamo is expected to be more efficient 
between the spiral arms ( Chamandy et al.|2015| l. 

The total magnetic held Btot can be presented as the quadratic sum of the reg¬ 
ular held Breg generated by the a — Q dynamo and the component Bturb generated 
by the small-scale dynamo. Tangling of Breg can produce another component, Btan- 
Btan and Bturb should be comparable, because both components are driven by the 
same turbulent gas motions. On the solar surface, the small-scale dynamo seems to 
be weak ( Stenho|20121 l, but may still operate in the solar interior. In contrast to the 
Sun, stars and planets, the regime of dynamo action is fully accessible in galaxies. 
The existence of the small-scale dynamo is of fundamental importance for the evolu¬ 
tion of magnetic helds in galaxies. Without a small-scale dynamo the seed helds in 
protogalaxies would be much weaker and the generation of large-scale regular helds 
would take much longer. Turbulent and tangled helds can be distinguished with (fu¬ 
ture) high-resolution observations. 


3 Observational tools 

Magnetic helds need illumination to be detectable. Polarized emission at optical, in¬ 
frared, submillimeter and radio wavelengths holds the clue to measure magnetic helds 
in galaxies. Huge progress has been achieved in the last decade. 

The various methods to measure magnetic helds are briehy discussed in the fol¬ 
lowing. The basic concepts are presented in more detail in [Klein & Fletcher] ( |2015| l. 

* Magnetic helicity H is defined as H = A- (V xA), where A is the magnetic vector potential. 
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3.1 Basic magnetic field components 

Table [^summarizes the different field components and the methods used to observe 
them. 


Table 1 Magnetic field components and methods to observe them. The symbols _L and || indicate field 
components perpendicular and parallel to the line of sight, respectively. 


Field component 

Notation 

Observational signatures 

Total field perpendicular 

^tot,± = ^turb.X +®reE.± 

Total synchrotron intensity 

to the line of sight 

Turbulent or tangled field 
perpendicular to the line of sight 

p2 p2 _|_ p2 

■°turb.± “ -^iso.i “‘""aniso.l 

Total synchrotron emission, 
partly polarized 

Isotropic turbulent or tangled 
field perpendicular 
to the line of sight 

.Siso.x {= \/2/3Biso) 

Unpolarized synchr. intensity, 
beam depolarization, 

Faraday depolarization 

Isotropic turbulent or tangled 

^iso.ll (— V^l/3-Siso) 

Faraday depolarization 

field along the line of sight 

Ordered field perpendicular 
to the line of sight 

p2 p2 1 p2 

■°ord.± “ "anisoa ^ "reg,^ 

Intensity and vectors of radio, 
optical, IR & submm pol. 

Anisotropic turbulent or 
tangled field perpendicular 
to the line of sight 

•^aniso._L 

Intensity and vectors of radio, 
optical, IR & submm pol., 
Faraday depolarization 

Regular field perpendicular 
to the line of sight 

^reg,_L 

Intensity and vectors of radio, 
optical, IR & submm pol., 
Goldreich-Kylafis effect 

Regular field along the 
line of sight 


Faraday rotation + depol., 
longitudinal Zeeman effect 


The total magnetic held is separated into a regular and a turbulent component. A 
regular held has a well-dehned direction within the telescope beam, while a turbulent 
held frequently reverses its direction within the telescope beam. In other words, the 
coherence scale of regular helds is much larger than the turbulent scale. Turbulent 
helds can be isotropic turbulent (i.e. the same dispersion in all three spatial dimen¬ 
sions) or anisotropic turbulent (i.e. different dispersions). Anisotropic turbulent helds 
originate from isotropic turbulent helds by the action of compressing or shearing gas 
Hows. 

The power spectra of turbulent helds are generally assumed to be of Kolmogorov 
type, with the largest scale given by the scale driving turbulence (e.g. supernova rem¬ 
nants in galaxies) and the smallest scale given by magnetic diffusivity. The power 
spectra of regular helds generated by the mean-held dynamo show a broad maximum 
over a range of scales which develops with time ( Moss et al.|2012| l. 

Present-day radio polarization observations with limited spatial resolution cannot 
distinguish turbulent helds generated by turbulent gas hows from tangled helds gen¬ 
erated from regular helds by small-scale gas motions; the components of both types 
of helds in the sky plane give rise to unpolarized synchrotron emission. 

Anisotropic turbulent, anisotropic tangled and regular held components contribute 
to the ordered held, observable in polarized synchrotron emission. Polarization angles 
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are ambiguous by 180°, so that polarized emission is not sensitive to field reversals. 
Faraday rotation and the longitudinal Zeeman effect are sensitive to the field direction 
and hence trace regular fields. 


3.2 Optical, infrared and submm polarization 


Elongated dust grains can be oriented with their major axis perpendicular to the inter¬ 
stellar magnetic field lines by paramagnetic alignment (the Davis-Greenstein effect) 
or, more efficiently, by radiative torque alignment ( Hoang & Lazarian||2008 2014| l. 
When particles on the line of sight to a star are oriented with their major axis perpen¬ 
dicular to the line of sight (and the field is oriented in the same plane), the different 
levels of extinction along the major and the minor axis leads to polarization of the 
starlight, with the E-vectors pointing parallel to the component of the interstellar 
magnetic field perpendicular to the line of sight (Eig.[2l. The interpretation of optical 
and near-infrared polarization observations of individual stars or diffuse starlight is 
based on this mechanism. Extinction is most efficient for grains of sizes similar to the 
wavelength. These small particles are only aligned in the medium between molecular 
clouds, not in the dense clouds themselves ( |Cho & Lazarian|20051 l. 

The detailed physics of the alignment is complicated and not fully understood. 
The degree of polarization p generated in a volume element along the line of sight 
of extent L is proportional to L), but also depends on the magnetic properties, 

density and temperature of the grains ( Ellis & Axon|1978) l. 

Starlight polarization yields the orientation of large-scale magnetic fields in the 
Milky Way (Sect. |5]l, and in the galaxies M 51 (jScarrott et al.||1987|l, NGC 6946 
(Fig .[T]) and the Small Magellanic Cloud ( |Lobo Gomes et al.|2015| l. The major short¬ 
coming when applying this method to extended sources like gas clouds or galaxies 
is that light can also be polarized by scattering, a contamination that is unrelated to 
magnetic fields and difficult to subtract. 

At far-infrared and submillimeter wavelengths, elongated dust grains emit lin¬ 
early polarized waves, without contributions by polarized scattered light. The B- 
vector is parallel to the magnetic field. The field structure can be mapped in gas 


clouds of the Milky Way, e.g. in the massive star formation site W 51 (Tang et al. 


20091 and in galaxies, e.g. in the halo of the starburst galaxy M 82 (Greaves et al. 
20001 . 


Magnetic force (total field strength B) and kinetic force by turbulent gas motions 
(turbulent velocity Vturb and density p) are in competition in the interstellar medium. 
In a strong field the field lines are straight and the dispersion of polarization angles CJ 
is small. According to Chandrasekhar & Fermi ( 1953|l, the total field strength Btot.x 
in the sky plane can be estimated from: 




( 2 ) 


Application to starlight polarization in the Milky Way yielded field strengths of about 
7 pG, in agreement with other methods (Sect.|^. The Chandrasekhar-Fermi method 
was improved by correcting for the observational errors in polarization angle and 






























Magnetic Fields in Spiral Galaxies 


7 


NGC6946 



Fig. 1 Optical emission (contours) and polarization E-vectors of NGC 6946, observed with the polarime- 
ter of the Landessternwarte Heidelberg at the MPIA Calar Alto 1-m telescope. The optical E-vectors are 
along the spiral arms in the eastern and western regions (compare with the radio polarization B-vectors in 
Fig.|13|. Polarization due to light scattering is probably occuning in the southern part (from|Fendt et al.| 
jl998t ) 


signal integration ( Hildebrand et al.|[2009 Houde etar]|2009 201 l| l. It can also be 
applied to radio polarization angles, like in M 51 (Houde et al.|2013li. 


3.3 Zeeman effect 


The Zeeman effect is the most direct method of remote sensing of magnetic fields. In 
the presence of a magnetic field B|| along the line of sight, a spectral line is split into 
two circularly polarized components of opposite sign (longitudinal Zeeman effect). 
The frequency shift is 2.8 MHz/Gauss for the Hi line and larger for molecular lines 
like OH, CN or H 2 O ( |Heiles & Crutche7||2005| l. The Zeeman effect is sensitive to 
regular fields, but turbulent fields can also be measured from the dispersion of the 


circularly polarized signal at a large number of locations on the source (Watson et al. 


20011. The interpretation of Zeeman splitting is hampered by its weakness, requiring 


high line intensities and careful correction of instrumental polarization. Most results 
have been obtained for the Hi line that traces clouds of diffuse (warm) gas in the 
Milky Way (|Crutcher et al.|2010|l and for the OH line for starburst galaxies (|Robishaw 
[etakp^. 


Magnetic fields perpendicular to the line of sight generate two shifted lines in 
addition to the main unshifted line, all linearly polarized (transversal Zeeman ef¬ 
fect). These lines cannot be resolved for observations in the Milky Way and external 
galaxies and no net polarization is observed under symmetric conditions. Detection 
of linearly polarized lines becomes possible for unequal populations of the differ- 
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ent sublevels, a gradient in optical depth or velocity, or an anisotropic velocity field 
( Goldreich & Kylafis||198l| l. Depending on the line ratios, the orientation of linear 
polarization can be parallel or perpendicular to the magnetic field orientation. This 
Goldreich-Kylafis effect was detected in molecular clouds, star-forming regions, out¬ 
flows of young stellar systems and supernova remnants of the Milky Way. It was also 
observed in the ISM of the galaxy M 33 ( Li & Henning|20rT] l, where it is consistent 
with the field orientations along the spiral arms as measured by polarized radio con¬ 
tinuum emission ( Tabatabaei et al.|[2008| l. The ±90° ambiguity needs to be solved 
with help of additional continuum polarization data. 


3.4 Total radio synchrotron emission and the equipartition assumption 


Radio continuum emission from galaxies is a mixture of thermal and nonthermal 
(synchrotron) components. The average thermal fractions are typically 10% and 25% 
at 20 cm and 6 cm wavelengths, respectively. At long wavelengths thermal emission 
is negligible, but thermal absorption may reduce the synchrotron intensity in star¬ 
forming regions. A first-order separation of the thermal and nonthermal components 
can be performed by assuming a uniform spectral index of a = 0.8 — 1.0 for the 
nonthermal intensity (/ v^“) and using the spectral index of the thermal emission 

in an optically thin plasma of a =0.1. However, synchrotron intensity depends on 
energy and age of cosmic-ray electrons (CREs) propagating from the spiral arms and 
losing energy by synchrotron radiation and inverse Compton scattering with photons, 
so that the synchrotron spectrum is expected to be flatter (smaller a) in spiral arms 
than in interarm regions, outer disk and halo. The advanced method to separate the 


two components introduced by Tabatabaei et al. (2007 1 uses an extinction-corrected 
Ha image as a template for the thermal radio emission and was first applied to the 
galaxy M 33. Simplified separation methods use an Ha image ( Heesen et al.|2014 l 
or a combination of Ha and 24 /rm infrared images as a thermal template that is 
subtracted from the total emission (|Basu & Roy|2013Tl. 


The intensity of total synchrotron emission (examples in Figs. 11 15 16 18- 


20 and 23 21^ all of which show mostly synchrotron emission) is a measure of the 
number density of cosmic-ray electrons (CREs) in the relevant energy range and of 
the strength of the total magnetic field component Btot,± perpendicular to the line of 
sight (i.e. in the sky plane). Synchrotron emission at a frequency v is emitted by a 
range of CREs with average energy E (in GeV) in a field of strength Btot,± (in MG)' 


V ~ 16MHz£^B, 


tot._L • 


(3) 


Equipartition is expected between the energy densities of the total cosmic rays 
(dominated by protons) and the total magnetic field. The cosmic-ray energy density is 
determined by integrating over their energy spectrum. This allows us to calculate the 
total magnetic field strength from the synchrotron intensity ( Beck & Krause|[2005 
Arbutina et al.|[20T^. The synchrotron intensity Lyn scales with the total magnetic 


field strength as Lyn Vice versa, Btot.x scales as: 


tot,_L 


B 


tot,± ' 


{hyn{Ko + l)/L) V(3+“), 


(4) 
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where a is the synchrotron spectral index (/gyn V^“) and L is the effective path- 
length through the source. Kq is the ratio of number densities of CR protons and 
electrons. Kq ~ 100 is a reasonable assumption in the star forming regions of the disk 
( Bell||1978l l. For an electron-positron plasma, e.g. in jets of radio galaxies, Kq — 0. 
The input parameters, especially the pathlength L and the ratio Kq are not well con¬ 
strained. Changing one of these values by a factor of e.g. 2 changes the field strength 
only weakly by2^/(^+“) ~ 1.2. 

Issues with the equipartition estimate are: 

(1) Eq.|^can be applied only for steep radio spectra with a > 0.5. For flatter spectra, 
the integration over the energy spectrum of the cosmic rays has to be restricted to a 
limited energy interval. 

(2) Due to the highly nonlinear dependence of /syn on Btot.x. the average equipar¬ 
tition value Btot.x derived from synchrotron intensity is biased towards high field 
strengths and hence is an overestimate if Btot varies along the line of sight or across 
the telescope beam ( Beck et al.|2003 1l. For a = \ and the equipartition case, the over¬ 
estimation factor g of the total field is (see Appendix A in Stepanov et al.|20f4|l: 


s = (<» 


>)‘'X<B»,,x> = (i + 5e^+5e‘)''X 
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tot,_L 


(5) 


where < > indicates the averages along the line of sight and the beam, and Q = 
(< ^ >yl'^ I < Btot,x > is the amplitude of the field fluctuations relative to the 

mean field. For strong fluctuations of 2 = 1, the overestimation factor is 1.46. 

(3) If energy losses of CREs are significant, especially in starburst regions or massive 
spiral arms, the equipartition values are lower limits and are underestimated. The ratio 
K increases because energy losses of aging CREs are much more severe than those 
of cosmic-ray protons. Using the standard value Kq underestimates the tot al magnetic 
held by a factor of {K/KqY !in such regions ( |Beck & Krause|2005| l. The same 
applies to the outer disk and halo of a galaxy, where the emitting CREs propagated 
far away from the places of origin and energy losses are significant. 

(4) In dense gas, e.g. in starburst regions, secondary positrons and electrons may 
be responsible for most of the radio emission via pion decay. Notably, the ratio of 
protons to secondary electrons is also about 100 for typical radio wavelengths (|Lacki 
|& Beck|201^ . 

(5) Energy equipartition needs time to develop and hence is not valid on small scales 
in time and space. A correlation analysis of radio continuum images from the Milky 
Way and the nearby galaxy M 33 by Stepanov et al. ( |2014| l showed that equipartition 
does not hold on scales smaller than about 1 kpc, which is understandable in view 
of the typical propagation length of CREs away from their sources (Sect. 4.3 i. The 
equipartition held strength is underestimated in regions around the CRE sources and 
overestimated in regions far away from the sources. 

Arguments for the validity of equipartition on large (kpc) scales come from the 
joint analysis of radio continuum and y-ray data, allowing an independent determina¬ 


tion of magnetic held strengths, e.g. in the Milky Way (Strong et al. 


20001, the Large 


Magellanic Cloud (EMC) ( |Mao et al.|[20T2] l and in M 82 ( |Yoast-Hull et al.|[20T3l ) 
Eurthermore, the estimates of strengths of large-scale ordered fields derived from 


but probably not in NGC 253 Yoast-Hull et al. 


2014 
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Faraday rotation data are similar to those from equipartition in several galaxies (see 
Table 2 in |Van Eck et al.|2015| l. 


3.5 Polarized radio synchrotron emission 


Linearly polarized synchrotron emission (examples in Figs. 0 13 [min] and 0 
emerges from CREs in ordered fields in the sky plane. As polarization angles are 
ambiguous by 180°, they cannot distinguish regular fields, defined to have a constant 
direction within the telescope beam, from anisotropic turbulent or tangled fields. 
Unpolarized synchrotron emission indicates isotropic turbulent or tangled fields that 
have random directions in 3D and have been amplified by turbulent gas flows (see 
Sect. [XT]). 

The intrinsic degree of linear polarization po of synchrotron emission in a per¬ 
fectly regular field is; 

Po = (1+ 0!)/(-+ a), (6) 

where a is the spectral index of the synchrotron emission. In spiral galaxies, typical 
values are a = 0.8 — 1.0, so that po = 73 — 75%. 

The observed degree of polarization is smaller than po due to the contribution 
from unpolarized thermal emission, wavelength-dependent Faraday depolarization 
(Sect. |3.7[ ) and wavelength-independent “beam depolarization” that is discussed in 
the following. 

Turbulent magnetic fields in galaxies preserve their direction over a coherence 
scale that depends on the properties of turbulence. 0 In the case of isotropic turbu¬ 
lent fields with a constant coherence length d, the observable degree of synchrotron 
polarization is (see |Sokoloff et al.] ( |1998| ) for a detailed discussion); 

p = poN-^f\ (7) 


where N ~ D^Lf/d^ is the number of synchrotron-emitting turbulent cells with di¬ 
ameter d and filling factor / within the volume defined by the telescope beam with 
the linear size D at the galaxy’s distance and pathlength L through the source. Eq. (|3 
is valid only at short radio wavelengths, where Earaday depolarization is small. 

Typical degrees of synchrotron polarization observed in nearby galaxies between 
3 cm and 6 cm wavelengths with typical spatial resolutions of D ~ 500pc are 1 — 5% 


in spiral arms (where isotropic turbulent fields dominate), e.g. in NGC 6946 (Beck 


2007[ ) and in M 33 ( Tabatabaei et al.|2008| l. With L ~ 1000 pc and / ~ 0.5 for the 


turbulent field coupled to the warm diffuse gas, we get d ~ 50 pc, consistent with 
estimates from Earaday depolarization at long wavelengths (Eq. 151. 

If the medium is pervaded by an isotropic turbulent held Bjso plus an ordered held 
Bord (regular and/or anisotropic turbulent) that has a constant orientation in the vol¬ 
ume observed by the telescope beam, the wavelength-independent polarization (for 


^ The degrees of polarization expected from compressed or sheared fields are discussed in the Appendix 
of|Beck et al.H2065|. 
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the case of equipartition between the energy densities of magnetic field and cosmic 
rays) amounts to (Sokoloff et al.|19^l: 


= pQ{\+3.5q^) / {\+A.5q^ +2.5q^) , 


P = po\ 

where q = Biso.x /Bord.x and Biso,x = \/^ 


( 8 ) 


3.6 Faraday rotation and Faraday Synthesis 


The polarization plane is rotated in a magnetized thermal plasma by Faraday rotation. 
As the rotation angle is sensitive to the sign of the field direction, only regular fields 
give rise to Faraday rotation, while the Faraday rotation contributions from turbulent 
fields largely cancel along the line of sight. Measurements of Faraday rotation from 


multi-wavelength observations (examples in Figs. 12 andp^ yield the strength and 
direction of the average regular field component along the line of sight. If Faraday 
rotation is low (in galaxies typically at wavelengths shorter than a few centimeters), 
the B-vector of polarized emission gives the intrinsic field orientation in the sky 
plane, so that the magnetic pattern can be mapped directly (Sect. [43| l. 

The Faraday rotation measure RM is defined as the slope of the observed variation 
of the polarization rotation angle % with the square of wavelength A; 


X = Xo + RMX^ 


(9) 


where Xo is the intrinsic polarization angle. A is measured in meters and RM in 
radm^^. RM is constant if is a linear function of A^, e.g. if one or more Faraday- 
rotating regions are located in front of the emitting region (Faraday screen). 

A nonlinear variation of % with A^ and hence a variation of RM with A^ occurs 
in case of: 

(1) emission and rotation in the same region where the distribution of electrons or 
regular magnetic field strength along the line of sight is not symmetric, or field rever¬ 
sals occur, or the field is helical ( [Sokoloff et al.|1998| l, 

(2) Faraday depolarization (Sect. 3.71, 


(3) several distinct emitting and rotating regions with different contributions to Fara¬ 
day rotation are located within the beam or along the line of sight or within the source 
(internal structure). 

In these cases, RM needs to be replaced by the Faraday depth 0 (Burn|1966l: 


observer 


A = Zo + 


0 A2 , where 0 = 0.812 J 


^r6gj| n^di , 


( 10 ) 


where 0 is measured in radm the line-of-sight magnetic field Breg,|| in MG’ the 
thermal electron density n^ in cm^^ and the line of sight I in pc. 

To measure the contribution of various sources to 0, wide-hand, multi-channel 
spectro-polarimetric radio data are needed. These are Fourier-transformed by a soft¬ 
ware tool called RM Synthesis ( jBrentjens & de Bruynj 200^ ^aM j^fS] ) or Faraday 
Synthesis (following the terminology proposed by Sun et al. 2015|l, to obtain the 
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(a) Sketch of different components along the line of sight. Some of them are 
emitting {Ei), are Faraday-rotating {Ri), or both. 
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(b) Resulting Faraday spectrum F{0), where 0 is the Faraday depth, of the 
components depicted in|2(a)| 


Fig. 2 Example for a set of different polarization-emitting and rotating components along the line of sight 
and the resulting Faraday spectrum F(0). El and £3 are point sources, E2,R2 is an extended emitting and 
rotating region. Rl, R3, and R4 are Faraday screens that are not emitting. Faraday depolarization effects 
are neglected here for simplicity (fromlGieBiihel et al.H2013)). 


Faraday spectrum F{<1>) (previously called Faraday dispersion function) at each po¬ 
sition of the radio image. F{^) shows the intensity of polarized emission and its 
polarization angle as a function of <l> (Fig.|^. 

If the medium has a relatively simple structure, Faraday spectra at several po¬ 
sitions can reveal the 3D structure of the magnetized interstellar medium (Faraday 
tomography). For example, regular fields, field reversals and turbulent fields can be 


recognized from their signatures in Faraday spectra (Bell et al. 2011 Frick et al. 


201 ![ Beck et al.|2012b|l. Helical fields can also imprint characteristic features in the 


Faraday spectrum ( Brandenburg & Stepanov|2014 Horellou & Fletcher|2014| l. Re¬ 

alistic galaxies have complicated Faraday spectra, as revealed from simplified galaxy 
models pdeguchi et al.|20r4| . 

Crucial to successful application of Faraday Synthesis is that the observations 
cover a wide range in A^, giving high resolution in Faraday spectra. Even for a wide 
frequency coverage, two components with a difference in intrinsic polarization angle 
of about 90° cannot be properly recovered by Faraday Synthesis, so that model fitting 
to the data in Stokes Q and U is needed ( [Farnsworth et al.|2011[ ). 

A grid ofRM measurements of polarized background sources is another powerful 
tool to measure magnetic field patterns in galaxies ( [Stepanov et al.|2008[ [Mao et ah] 
[2012[). At least 10 RM values from sources behind a galaxy’s disk are needed to 
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recognize a simple large-scale field pattern if the Galactic foreground contribution is 
constant and the background sources have no intrinsic contributions to RM. 

Measuring the strength of the regular field from RM needs additional information 
about the thermal electron density, e.g. from the arrival times of pulsar signals that 
are delayed in a cloud of ionized gas proportional to the dispersion measure DM = 
f n^dl. The ratio 


RM/DM = 0.812 <B,eg ||ne > / < He > (H) 

allows us to compute the field strength Breg.|| (in MG) if the fluctuations in field 
strength and in electron density are uncorrelated. The value for Breg,|| derived from 
pulsar data is an underestimate if small-scale fluctuations in field strength and in 
electron density are anticorrelated, as expected for local pressure equilibrium, while 
it is an overestimate if the fluctations are correlated, as expected for compression by 
shock fronts or turbulence (|Beck et al.|2003|l. 


3.7 Faraday depolarization 


In a single region containing CREs, thermal electrons and purely regular magnetic 
fields, wavelength-dependent Faraday depolarization occurs because the polarization 
planes of waves from the far side of the emitting layer are more rotated than those 
from the near side. This effect is called differential Faraday rotation. For one single 
layer with a symmetric distribution of thermal electron density and field strength 
along the line of sight the degree of polarization is reduced to (|Burn|1966[ |Sokoloff 
|etal.|1998D : 


p = PQ \sin{2RMX^)/{2RMX^)\, 


( 12 ) 


where RM is the observed rotation measure from the source, which is half of the 
Faraday depth <l> through the whole layer, p varies periodically with the square of 
wavelength. 


Applying Faraday Synthesis (Sect. 3.6 1 to polarization data of a region that is 


emitting and Faraday-rotating reveals an extended component in the Faraday spec¬ 


trum (Fig. 1^. However, regions broader than ^>max — TtjX^ 
where Amin is the minimum wavelength of the observations ([Brentjens & de Bruyn 


cannot be recovered. 


|2005| ); only two “horns” remain visible at the edges of the structure in the Faraday 
spectrum F{<P) (Beck et al. 2012bI. This problem is similar to the missing short 
baselines in synthesis imaging. 

For multiple emitting + rotating layers, the Faraday spectrum contains several 
extended components and Eq. ( [T^ is no longer applicable. Shneider et al. (20141 
extended Burn’s model to a multi-layer (e.g. galaxy disk + halo) medium. 

Faraday rotation in helical fields has a completely different behaviour and may 
lead to an increase of the degree of polarization with increasing wavelength (“re¬ 
polarization”) in certain wavelength ranges ( |Sokoloff et al.|[T998[ [Brandenburg &| 
|Stepanov|2014||Horellou & Fletcher|2014| ). 
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Turbulent fields also cause wavelength-dependent depolarization, called Fara¬ 
day dispersion ( |Sokoloff et al.|19^ Arshakian & Beck|2011| l. For an emitting and 
Faraday-rotating region (internal dispersion) the degree of polarization is reduced to: 


F = Fo (1-exp(-5))/5, (13) 

where S = Orm is the dispersion in rotation measure and depends on the 

turbulent field strength along the line of sight, the turbulence scale, the thermal elec¬ 
tron density, and the pathlength through the medium. 

Depolarization by external dispersion occurs in a turbulent Faraday-rotating (but 
not emitting) medium in the foreground if the diameter of the telescope beam at the 
distance of the screen is larger than the turbulence scale: 

p = poexp{-S). (14) 


An important effect of Faraday dispersion is that the interstellar medium becomes 
“Faraday thick” for polarized radio emission beyond a certain wavelength, depending 
on Orm, and only a front layer remains visible in polarized intensity. Typical values 
for galaxy disks are Orm = 20 — 100 radm~^, while star-form ing regions can have 
dispersions of cJrm — 1000radm^^ ( Arshakian & Beck|201 11 . The maximum in the 
spectrum of polarized intensity falls into the wavelength range 2-20 cm. 

1/2 

In a random-walk approach, we may write Orm = 0.812Biso || ngdN^^ (Beck 
20071), where Bisn ii = y/l/S Fiso is the strength of the isotropic turbulent field, n 


the electron density within a cell of size d, and A|| = Lf /d is the number of cells 
along the line of sight L with a volume filling factor /. The average electron density 
along the line of sight is < Ue >= «e//, so that we get: 


Crm = 0 . 812 Bis 


< He > 


^/Ldff. 


(15) 


The average value for galaxy disks of Orm ~ 50 rad m 
values for the warm diffuse ISM of < n^ >~ 0.03 cm^^ 
Berkhuijsen et al.|2006 1 , Fiso — lOpG, L ~ 1000pc, and d ~ 50pc. 


^ is consistent with typical 
and / ~ 0.5 (Fig. 10 in 


The dispersion Orm leads to a “Faraday forest” of N components in the Faraday 
spectrum. If N is not large, the components are possibly resolvable with very high 


Faraday resolution, hence a wide span of the observations (Beck et al. 2012b 
|Bemet et al.|2012| ). 


4 Results 

4.1 Magnetic field strengths 

Strength is the fundamental quantity to estimate the dynamical importance of mag¬ 
netic fields. The equipartition assumption (Sect. [3A| is the only presently applicable 
method to measure the total field strength. Zeeman measurements will become possi¬ 
ble in spiral galaxies with future radio telescopes, but will trace fields on small scales 
only. 
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The average equipartition strength (Eq.|^ of total fields (corrected for inclination) 
for a sample of 74 spiral galaxies is Btot = 9±2/rG (Niklas 19951. The average 
strength of 21 bright galaxies observed between 2000 and 2010 is Btot = 17 ± 3 /iG 
( Fletcher||2010| l. Gas-rich spiral galaxies with high star-formation rates, like M 51 
(Fig. 1^, M 83 (Fig. 1^ and NGC 6946 (Fig. [T3] l, have total field strengths of 20 — 
30/iG in their spiral arms. The strongest total fields of 50— 100/iG are found in 
starburst galaxies like M 82 (|Adebahr et al.||2013|), in the “Antennae” NGC 4038/9 
(|Chyzy & Beck|[2004|l, in nuclear starburst regions like in NGC 253 (Heesen et al. 


201 la I, and in barred galaxies ([Beck et al.|2005|l. In such galaxies, energy losses of 


CRFs can be strong, so that the equipartition values are underestimates (Sect. 3.4 1 . 
Radio-faint galaxies like M 31 (Fig. [nj and M 33 have weaker total magnetic fields 
(about 6 /rG). The similarity to the values derived for our Milky Way with various 
methods (Sect. gives confidence that the equipartition assumption is valid. 

Field strengths of 0.5 — 18mG were detected in starburst galaxies by the Zee¬ 


man effect in the OH megamaser emission line at 18 cm wavelength (Robishaw et al. 


20081. These values refer to highly compressed gas clouds and are not typical of the 


diffuse interstellar medium. 

The strength of ordered fields in the sky plane as observed by polarized syn¬ 
chrotron emission varies strongly between galaxies, from 10-15 /iG in M 51 and the 
magnetic arms of NGC 6946 (Fig. [T3]l to about 5 /iG in the star-forming ring of M 31 
(Fig.[TT]i. The average strength of the ordered fields of 21 bright galaxies observed 
since 2000 is Bord,± = 5 /iG with a standard deviation of 3 /iG and an average ratio 
/^ord.x//^tot of 0.3 (Fletcher 20101. 

The only deep observation of an Sa galaxy, M 104, with a prominent dust ring, 
revealed ordered, partly regular magnetic fields ( Krause et al.|2006] l. Spiral galaxies 
of type SO and elliptical galaxies without an active nucleus have very little star forma¬ 
tion and hence hardly produce CRFs that could emit synchrotron emission. Irregular 
starburst galaxies show only spots of ordered fields ( |Heesen et al.|201 lb) !. No ordered 
fields could be detected so far in dwarf irregular galaxies (Chyzy et al.|201 l]l. 


4.2 Magnetic energy densities 


The relative importance of various competing forces in the interstellar medium can be 
estimated by comparing the corresponding energy densities. The mean energy densi¬ 
ties of the total (mostly turbulent) magnetic field and the total cos mic rays, av eraged 
in rings of about 1 kpc width, are ~ 10^ " erg cm^^ in NGC 6946 (Beck 2007 i, M 63, 
M 83, NGC 4736 (|Basu & Roy|2013| and IC 342 (Fig.[|), and ~ 10^^^ erg cm^^ in 
M 33 ( Tabatabaei et al.|2008 i. In all these galaxies the average magnetic energy den¬ 
sity is similar to the average kinetic energy density of the turbulent motions of the 
neutral gas across the star-forming disk (assuming a constant turbulent velocity of 
lOkms^*) and about 10 times larger than that of the warm ionized gas (but still 500- 
1000 times smaller than the energy density of the general rotation of the neutral gas). 
This result is similar to that derived for the Milky Way ( |Cox|2005] l. The ISM is a low- 
P plasma (where j3 is the ratio of thermal to magnetic energy densities). Magnetic 
fields are dynamically important. 
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Fig. 3 Radial variation of the energy densities in IC 342, determined from obsei'vations of synchrotron and 
thermal radio continuum and the CO and Hi line emissions: for the total magnetic field (black tri¬ 
angles), identical to that for the total cosmic rays, for the ordered magnetic field (green triangles), 

for the turbulent motion of the neutral gas £turb = 0.5 p„ v^rb^ (blue squares), where Vturb = 2 km/s, and 
for the thermal energy of the warm ionized gas fta, = 1.5 < Me > (red circles), where Te = 10"* K. The 
en'or bars include only errors due to rms noise in the images from which the energy densities are derived. 
No systematic errors are included, e.g. those imposed by a radial variation of thermal gas temperature, 
filling factor or turbulent gas velocity, nor enors due to deviations from the equipartition assumption (from 

[BeckUM^ ). 


According to Fig.[^ Fturb/Fth ^ F i-e. turbulence is supersonic. Supersonic tur¬ 
bulence leads to shocks and hence dissipation of kinetic energy into heat, until turbu¬ 
lence becomes transonic and the system reaches an equilibrium state. The observa¬ 
tion of supersonic turbulence may indicate that a quasi-steady state is never reached. 
However, observational bias effects need to be considered. The energy density of hot 
gas in the ISM, neglected in Fig. is similar or somewhat larger than that of warm 
gas, depending on its volume filling factor (Ferriere 200 l| l, so that its inclusion would 
not change the above result significantly. Another bias is the assumption of a constant 
gas temperature T^. If cosmic rays propagate via the streaming instability (Sect. 4.3 1 , 
they may heat the gas and increase in the outer disk ([Wiener et al.|2013]l. 


Another important message from Fig. (and similarly from the other galaxies 
studied so far) is that the magnetic energy density decreases radially more slowly 
than the kinetic energy density and dominates in the outer galaxy disk, i.e. the Alfven 
velocity va Btot/ \fp is larger than the turbulent velocity and increases with ra¬ 
dius, which may cause magnetic effects on the rotation curve of the gas (jElstner et al. 


outer galaxy disk than shown in Fig. because the equipartition field strengths are 
underestimates due to energy losses of the CREs (Sect. [374)1. Eurthermore, the turbu- 


20141. Note that the ratio of magnetic to kinetic energy may be even higher in the 
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lent velocity tends to decrease with radius ( |Tamburro et al.|2009) l, which would also 
enhance the ratio of magnetic to kinetic energy. 

Possible reasons why equipaitition between magnetic and kinetic energy density 
does not hold within galaxies are: 

(1) energy or pressure balance is valid between the magnetic field and the sum of 
pressures related to all gas components (kinetic and thermal) ( Basu & Roy|2013] l; 

(2) the efficiency of the small-scale dynamo increases from the inner to the outer 
galaxy; 

(3) the magneto-rotational instability (MRI) leads to super-equipartition fields (Fig. 7 
in |Gressel et al.|2013| l; 

(4) the equipartition field strengths are underestimated in the inner galaxy and over¬ 
estimated in the outer galaxy due to CRE propagation (see Sect. |3.4| l, hence flattening 
the profile of the magnetic energy density. 

Scenario (4) is supported by the analysis of radio and IR images (|Berkhuijsen 


|et al.|2013l l and the cross-correlation based on wavelet transforms that breaks down 
below a certain scale that can be interpreted as the propagation length of CREs 
(|Tabatabaei et al.|2013a|l. CRE propagation in IC 342 seems to be exceptionally fast 
(Sect.|4.3». 

The energy density of the ordered magnetic field (Eig. even increases radially. 
This field may be assumed to be mostly regular, so that we can invoke the a — Q 
dynamo that also operates in the outer disk of galaxies ( Mikhailov et al.|2014 l. Al¬ 
though the star-formation activity is low in the outer disk, the magneto-rotational 
instability (MRI) may serve as the source of turbulence required for a — Q dynamo 
action ( Sellwood & Balbus]|1999 Gressel et al.||2013 l. Alternatively, regular fields 
generated in the inner disk could be transported outwards by the joint action of a 
dynamo and turbulent diffusivity (Moss et al.|1998|l. 


4.3 Cosmic-ray propagation 


The measurable extent Lsyn of synchrotron disks or halos (see Sect. 4.13 i gives us 


information about the extent of the total magnetic field and the propagation of CREs. 
Under energy equipartition conditions, the radial exponential scale length of the total 
field in a disk of mildly inclined galaxies (or the vertical exponential scale height 
in a halo of almost edge-on galaxies) is Lb ^ (3 + CtjLsyn. where a ~ 0.9 is the 
synchrotron spectral index. The scale length Lcr (or scale height) of the total cosmic 
rays is half of this value. These are lower limits because the CREs lose their energy 
with distance from the star-forming disk, so that the equipartition assumption yields 
too small values for the field strength (Sect. |3.4| l. The scale lengths of synchrotron 
disks of typically — 3 — 5kpc at 20 cm wavelength (Beck||2007 Basu & Roy 
2013t Mulcahy et al. ^14| l yield Lb — 12— 19 kpc and Lcr — 6 — 9.5 kpc. The nearby 
galaxy IC 342 (Eig. ^ r eveals an exceptionally large synchrotron scale length of 
about 16 kpc ( Becl^20T5] l. 

The radial scale length and its frequency dependence give us a general idea about 
the speed and type of CRE propagation. At wavelengths of > 1 m, synchrotron and 
inverse Compton losses are weaker than at shorter wavelengths, so that CREs can 
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diffuse further outwards from their locations of origin. This is clearly seen in M 31 
and M 33 ( Berkhuijsen et al.|2013| l. The LOFAR image of the galaxy M 51 at 2 m 
wavelength reveals a steepening of the radial distribution of synchrotron emission, 
located at about lOkpc radius, beyond the sharp decrease in the star-formation rate, 
indicating a propagation length of a few kpc ( Mulcahy et al.|[2OT4 l. The frequency 
dependence of the scale length favours diffusive propagation in the cases of M 51 
and several other galaxies ( |Basu & Roy|2013]|Berkhuijsen et al.|2013[|Mulcahy et al.| 
20141 and fast convective transport in IC 342 ( Beck|2015]r 


The galactic fields may extend further out into intergalactic space than visible 
in radio synchrotron images. A large radial scale length means that magnetic fields 
may affect the global rotation of the gas in the outer parts of spiral galaxies, possibly 


explaining some fraction of the flattening of the rotation curves (Battaner & Florido 
[20U7l[Elstner et al.|2014D . 


4.4 Relations between gas, star formation and magnetic fields 


Radio continuum and infrared emissions from galaxies are closely correlated. 

Firstly, a global correlation exists between the galaxy-integrated luminosity of the 
total radio continuum emission at around 20 cm wavelengths (frequencies of around 
1 GHz), which is mostly of synchrotron origin, and the infrared (IR) luminosity of 
star-forming galaxies. This is one of the tightest correlations known in astronomy. 
The correlation extends over five orders of magnitude ( Bell| 2003]l, is slightly nonlin¬ 
ear in log-log scale ( Bell|2003 Pierini et al.|2003|l with an exponent of 1.09 ± 0.05 
( Basii^t^alJ20^ and is valid in starburst galaxies to redshifts of at least four (Sey- 


mour et al.|2008|l, without evidence for evolution with redshift up to 1.2 (|Basu et ~ 


20151. As the correlation between radio thermal and IR luminosity should be strictly 
linear, the exponent of the nonthermal (synchrotron)-IR correlation is steeper than 
1.1. A breakdown of synchrotron emission and of its correlation with IR is expected 
beyond a critical redshift when inverse Compton loss of the CREs dominates syn¬ 
chrotron loss; this critical redshift will give us information about the held evolution 
in young galaxies ( [Schleicher & Beck|2013| l. 

In galaxies with strong magnetic fields, e.g. in starburst galaxies, the CREs may 
lose their energy by synchrotron loss within the galaxy disk (“lepton calorimeter 
model”, Volk ( 1989| l), where the correlation is linear and due to the joint generation 
of CREs and dust-heating UV photons in star-forming regions, independent of the 
magnetic held. In Milky Way-type galaxies, magnetic fields are crucial to understand 
the correlation. [Niklas & Beck| ( |1997j l proposed an “equipartition model”, where total 
(mostly turbulent) magnetic fields, cosmic rays, gas density and star formation are re¬ 
lated, and predicted a nonlinear synchrotron-IR correlation with an exponent of about 


1.3, in good agreement with observations. The more detailed analysis by Lacki et al. 


(2010 | l identified multiple feedback mechanisms that are not yet fully understood. 

Secondly, the total radio and IR intensities within galaxies are also highly corre¬ 
lated, but with a smaller exponent, due to CRE propagation ( Berkhuijsen et al.|2013| l. 
CRE propagation is also responsible for the flattening of the profile of cosmic-ray 
energy density in Eig. The exponent of the correlation was found to be different in 
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Fig. 4 Con'elation between the strength of the total, turbulent and ordered helds and the star-formation rate 
per area (determined from the 24 (im infrared intensities) within the galaxy NGC 6946. The points for 5tot 
and 5ord ^6 shifted by +0.4 and -0.4, respectively, to fit into the same figure. The solid lines are bisector 
fits with slopes of 0.14 ± 0.01 for Btot and 0.16 ± 0.01 for 5tuj-b, the dashed lines ordinal^ least-squares fits 
with slopes 0.089 ± 0.004 and 0.100 ± 0.005 (from [Thbatabaei et al.H2013b^ and priv. comm.). 



the central region, spiral arms and interarm regions in several nearby galaxies (Dumas 
let aLpOlltfEasu et al.|20l^ , and it also differs between galaxies due to differences 
in the diffusion length of CREs ( Berkhuijsen et al.|2013[ Tabatabaei et al.|2013a| l. 

The magnetic field and its structure play an important role in understanding the 
correlation ( Tabatabaei et al.|[2013a|b| l. The synchrotron-lR correlation can be pre¬ 
sented as a correlation between total field strength Btot and star-formation rate sur¬ 
face density Esfr with an exponent of 0.18 ± 0.01 within NGC 4254 ( |Chyzy|2008| l, 
0.14 ±0.01 within NGC 69 46 (Fig. |^, and 0.3 0 ±0.02 (globally and locally) for a 
sample of 17 spiral galaxies ( |Heesen^ al. 2014 1 . A very similar exponent is expected 
from theoretical considerations ( |Schleicher & Beck|2013| l. 

The Kennicutt-Schmidt law ±sfr valid within galaxies, with A ~ 1.0 

for the molecular gas and A ~ 1.5 for the atomic gas (Bigiel et al. 2008 f p] Then the 
above relation Btot -^sfr*^ corresponds to Btot if molecular gas dominates, 

e.g. in the inner regions of most spiral galaxies, and Z?tot -^gas^^ '^he outer regions 
that are dominated by atomic gas. 

Turbulent fields in spiral arms are probably generated by turbulent gas motions 
related to star formation activity (Sect.j^. The small-scale dynamo predicts equiparti- 
tion between magnetic and kinetic energy densities (Btot P^'^, assuming Vturb=const), 


which is supported by the observed average values of a galaxy sample (|Nildas & Beck 


^ The nonlinearity of the Kennicutt-Schmidt law can be interpreted as a dependence of the star- 
formation efficiency (SFE) on the stellar mass surface density jShi et al.|20n) . Interestingly, effects of 
the magnetic field on the SFE have not been investigated so far. 
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|1997[). H owever, within a galaxy the ratio between these energy densities decreases 
(Fig. [^, which is consistent with an exponent of < 0.5 found from the Zsfr data 
(see above). Possible reasons why equipartition between magnetic and kinetic energy 
density holds on average between galaxies, but not within galaxies, were discussed 
in Sect. 14.21 

The correlation between the radio continuum luminosity and the molecular gas 
(traced by its CO line emission) is similarly tight as the radio-IR correlation and is 
also nonlinear, with a slope of 1.31 ±0.09 ( |Liu & Gao|2010| l. The correlation between 
the synchrotron radio continuum and CO intensities, observed with 60 pc resolution 
within the spiral galaxy M 51 ( Schinnerer et al.|2013| l, appears to be tighter than the 
radio-IR correlation and may be the fundamental one. Much of the molecular gas in 
the spiral arms of M 51 is not directly associated with ongoing massive star forma¬ 
tion, implying that star formation cannot be the sole cause for the radio-CO relation. 
Either the synchrotron emission emerges from secondary CREs that are produced in 
the interaction of cosmic-ray protons with the dense molecular material ([Murgia et al. 


20051, or the synchrotron emission is from primary CREs in strong magnetic fields 


coupled to the dense gas (Nrklas & Beck 19971, in connection with fast diffusion 
of CREs along the spiral arms. High-resolution observations of nearby galaxies (e.g. 
with SKA and ALMA) can provide the necessary information on the coupling be¬ 
tween magnetic fields and gas on small scales, to test which of the proposed origins 
for the correlation is indeed true. 

In contrast to the turbulent field, the ordered field is either uncorrelated with 
the star-formation rat e (Fig. or antic orrelated in interarm regions, where the star- 
formation rate is low (|Frick ^al.|2001|l. 


4.5 Spiral fields 


The most striking result from the observation of ordered (regular and/or anisotropic) 
fields with help of polarized synchrotron emission is that spiral patterns are found 


in almost every spiral galaxy (Beck & Wielebinski 20131, in galaxies with a star¬ 
forming ring ( Chyzy & Buta|2008| , in flocculent galaxies without gaseous spiral arms 
( Soida et al.|2002| l, in central regions of galaxies (Fig.p^, and in circum-nuclear gas 
rings of barred galaxies (Fig.fThll. 

In galaxies with strong density waves like M 51 (Fig. and M 83 (Fig. en¬ 
hanced ordered fields occur at the inner edges of the inner optical arms where the 
cold molecular gas is densest ( |Patrikeev et al.|2006| l. These are probably anisotropic 
turbulent fields generated by the compression of the density wave. An analysis of 
dispersions of the radio polarization angles at 6 cm wavelength in M 51 shows that 
the ratio of the correlation lengths parallel and perpendicular to the local ordered 
magnetic field is about two ( Houde et al.|2013) . 

In galaxies without strong density waves the ordered fields are not coinciding 
with the gaseous or the optical spiral arms. The typical degree of radio polarization at 
small wavelengths is only a few percent within the spiral arms, hence the field in the 
spiral arms must be mostly tangled or randomly oriented within the telescope beam 
on scales of a few 10 pc (Eq.|^. 
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Fig. 5 Polarized radio emission (contours) and B-vectors of M 51, combin ed from observations a t 6 cm 
wavelength with the VLA and Effelsberg telescopes at 8" resolution (from [Hetcher et ^ 201 1}) . The 
background colour image shows the CO line emission from molecular gas (from|Helfer et al. j2003)). 


Ordered fields are strongest (10—15 /rG) in the regions between the spiral arms 
and oriented parallel to the adjacent spiral arms. In some galaxies they form magnetic 
arms, like in NGC 6946 (Fig. [T3] l, with exceptionally high degrees of polarization up 
to 50%. These are probably generated by the a — Q dynamo (Sect. [4!9] ). 

The ordered magnetic field in the galaxy IC 342 reveals several spiral polarization 
arms of different origins (Fig.[^. In contrast to NGC 6946, there is only a rudimentary 
magnetic arm in an interarm region in the north-west, probably because of weaker 
dynamo action in IC 342 (Sect. 4.9 1 . A narrow polarization arm of about 300pc 
width, displaced inwards with respect to the inner arm east of the central region by 
about 200pc, indicates that magnetic fields are compressed by a density wave, like in 
M 51. A broad polarization arm of 300-500 pc width around the northern optical arm 
shows systematic variations in polarized emission, polarization angles and Faraday 
rotation measures on a scale of about 2kpc, indicative of a helically twisted flux tube 
generated by the Parker instability (Sect. 4. lOi. Several broad polarization arms in 
the outer galaxy are coincident with spiral arms in the total neutral gas. 

At wavelengths of around 20 cm, most of the polarized emission from the far side 
of the disk and halo is Faraday-depolarized and the emission from the front side dom¬ 
inates. A striking asymmetry of the polarized emission occurs along the major axis 
of 12 spiral galaxies with inclinations of less than about 60°. The emission is always 
much weaker around the kinematically receding side (positive radial velocities) of 
the major axis (Urbanik et al. 1997[ Braun et al. 2010[ Vollmer et akpoO l. This 
asymmetry is still visible at 11 cm wavelength, but disappears at smaller wavelengths 
where the polarized emission from the far side becomes observable. In strongly in¬ 
clined galaxies, both sides of the major axis become Faraday-depolarized at around 
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Fig. 6 Polarization B-vectors of IC 342, combined from obser vations at 6 c m wavelength with the VLA 
and Effelsberg telescopes and smoothed to 25" resolution (from [Be^j2015t ), overlaid on a colour image 
from the Kitt Peak Observatory (credit: T.A. Rector, University of Alaska Anchorage, and H. Schweiker, 
WIYN and NOAO/AURA/NSF). A region of 16' x 16' (about 16 x 16kpc) is shown. (Copyright: MPIfR 
Bonn) 


20 cm. Modeling of a combination of disk and halo fields, as predicted by a — 
dynamo models (Sect. 2), can explain the asymmetry ( Braun et al.|2010| . 

At even longer wavelengths, Faraday effects depolarize the synchrotron emission 
almost completely. With help of Faraday Synthesis applied to 90 cm data from the 
Westerbork Synthesis Radio Telescope, an extremely low average degree of polariza¬ 
tion of 0.21 ±0.05% was measured in the star-forming “ring” of M 31 ([GieBubel et al. 


20131. No polarized emission could be detected from M 51 with the Low Frequency 
Array (LOFAR) at around 2 m wavelength - total Faraday depolarization occurs at 
such long wavelengths (|Mulcahy et al.|20T4)l. 


4.6 Spiral pitch angles 

In mean-field dynamo theory, the average pitch angle of the spiral magnetic field pQ 
and the average pitch angle of the gaseous spiral arms p^ are not related, while they 
are observed to be roughly similar in a limited sample of galaxies (Fig. |^. Within 
galaxies, too, the magnetic pitch angles are similar to the pitch angles of the gaseous 
spiral arm structures and of those of the polarized intensity, tracing the structures 
of the ordered field (e.g. Fig. |^. In IC 342 the magnetic pitch angle |pb| decreases 
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M83 6.2cm VLAtEffelsberg Polarized Intensity ■» B-Vectors HPBWsIS" 




Fig. 7 Polarized radio emission (contours) and 
B-vectors of M 83, combined from observations 
at 6 cm wavelength with the VLA and Effelsberg 
telescopes and smoothed to 15'^ resolution (from 
|Frick et ^j2Q16) ), overlaid onto an optical im¬ 
age from Dave Malin, Anglo Australian Obser¬ 
vatory. (Copyright: MPIfR Bonn and AAO) 


Fig. 8 The power spectra of pitch angles of structures 
at the wavelet-filtered scale F (about 2.6 kpc) in M 83, 
normalized to unit maximum: the pitch angles of the 
structures of polarized intensity at 6 cm wavelength 
(black dotted), total neutral gas (green) and molecular 
gas CO (blue dashed). The spectrum of the magnetic 
pitch angles is shown with the red dot-dashed curve 
(from [Fnck et al■|f2Q16^ )■ 


with radius (Fig.[T0|), following the pitch angle of the gaseous spiral structures (|Beck 
[20T5l ). 

The formation of spiral arms and spiral magnetic patterns seems to be related, as 
expected e.g. for MHD density waves (Lou et al. 1999| l. Almost identical pitch angles 
would indicate field alignment with the gas structures by shear or compression in 
density waves. However, based on anisotropic wavelet transformations, a systematic 
shift was detected in the barred galaxy M 83 (Fig.|^, in the sense that |pb| is larger 
than |pa| by about 20° (Fig. |^. Similar results were found in other galaxies. Such 
differences can be regarded as a signature of a — dynamo action (Sect. |4.8|l. 


In the thin-disk approximation of the a — Q dynamo, the magnetic pit ch angle is 
given by |pb| = where Ra and Rq are the dynamo numbers (Shukurov 


20051. Simplified estimates for the dynamo numbers and a flat rotation curve, which 


is valid beyond a radius of about a few kpc in most spiral galaxies, give |/7 b| — d/H, 
where d is the turbulence scale (about 50 pc, see [Fletcher et al. p011| l), and H is 
the scale height of the ionized gas. The scale heights of Hi disks in spiral galaxies 
increase radially ( |Bagetakos et al.|2011[ ). One may speculate that the disk of ionized 
gas also flares (e.g. assumed in the model by Gressel et al.|2013] l, which would ex¬ 
plain the decrease of \pQ \ observed in many galaxies. If so, the radial variation of \p^ \ 
in IC 342 (Fig. [T0| indicates that the scale height H is approximately constant until 
12' (about 12 kpc) radius and then increases. Alternatively, the magnetic pitch angle 
may be affected by gas flows, for instance by outflows that become weaker towards 
the outer disk and can decrease the effective dynamo number ([Shukurov et al.|2006|l. 
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Fig. 9 Average magnetic pitch angles |pb| and av¬ 
erage pitch angles of gaseous spiral stnrctures |pa| 
in several radial regions of 12 galaxies. The dashed 
line corresponds to |pb| = I Pa I (fromlVan Eck et al.l 


T 



Fig. 10 Azimuthally averaged magnetic pitch an¬ 
gle pb as a function of radius (in arcminutes) in the 
plane of IC 342. pB is negative in this galaxy, be¬ 
cause the spiral pattern (se e Fig.[^ turn s outward in 
the clockwise sense (from[Beck|j|2015^). 


4.7 Regular fields 


Spiral fields can be generated by compression in spiral arms, by shear in interarm re¬ 
gions or by dynamo action (Sect.|^. Measuring Faraday rotation (Sect. 3.6 1 is crucial 
to distinguish between these mechanisms. Large-scale sinusoidal patterns of Faraday 
rotation measures (RM) along azimuthal angle in the galaxy plane (modes) are sig¬ 
natures of regular fields generated by the a~Q dynamo and can be identihed in RM 


derived from polarized emission from the galaxy disks at several frequencies (Krause 


19901 or in RM data of polarized background sources ( [Stepanov et al.|2008) l. If sev¬ 
eral modes are superimposed, a Fourier analysis of the RM variation with azimuthal 
angle is needed. The resolution and sensitivity of present-day radio observations is 
sufficient to identify 2-3 modes. 

The disks of about a dozen nearby spiral galaxies reveal large-scale RM patterns. 
The Andromeda galaxy M 31 (Fig. [n]) is the prototype of a dynamo-generated ax- 
isymmetric spiral disk held, with a striking sinusoidal RM variation along the star¬ 
forming ring (Fig. JSL which is a strong indication of an axisymmetric spiral held 
(mode m = 0) ( [Fletcher et al.|2004 i. Other candidates for a dominating axisymmet¬ 
ric disk held are the nearby spirals IC 342 (Sokol off et al.[[T992[ Beck[[2015] l and 
NGC 253 ( Heesen et al.[2009b| l. The axisymmetric held in the irregular Large Magel¬ 
lanic Cloud (LMC) is almost azimuthal (i.e. small pitch angles) ( [Gaensler et al.[2005 ] 
Mao et al.[2012[ ). Dominating bisymmetric spiral helds (mode m = 1) are rare, as pre¬ 
dicted by dynamo models, but possibly exist in M 81 ([Krause et al.|1989[ [Sokoloff 
et al.|1992[ i. Faraday rotation in NGC 6946 and in other similar galaxies with mag¬ 
netic arms can be described by a superposition of two azimuthal modes (m = 0 and 
m = 2) with about equal amplitudes, where the quadrisymmetric (m — 2) mode is 
phase shifted with respect to the density wave ( [Beck|2007] l. For several other galax¬ 
ies, three modes (m = 0, 1 and 2) are necessary to describe the data (Table|^. 

The spiral pattern of magnetic helds cannot be solely the result of a —Q dynamo 
action. If the beautiful spiral pattern of M 51 seen in radio polarization (Fig. were 
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Table 2 Decomposition of regular magnetic fields in galaxies into azimuthal modes of order m. Columns 
2 — 4 give the approximate amplitudes relative to the dominating mode. A dash indicates modes that were 
not sought, whereas a zero indicates that no signature of this mode was found. Modes of order higher than 
2 cannot be detected with present-day telescopes. This is an updated version of a table in|Fletch^j2Q10^. 


Galaxy 

m — 0 

m—l 

m — l 

Reference 






M31 

1 

0 

0 

Fletcher et a] 

20041 




M33 

1 

1 

0.5 

I'abatabaei et 

al. 

2Dh8 




M51 (disk) 

1 

0 

0.5 

Fletcher et al 

2l 

m 




M51 (halo) 

0 

1 

0 

Fletcher et al 

■ 2Dm 




M81 

0.5 

1 

- 

iiolcolott et a 

1992 




NGC 253 

1 

- 

- 

Fleesen et al. 

Jll009b 




NGC 1097 

1 

1 

1 

Beck et al. 2 

0D5 





NGC 1365 

1 

1 

1 

Beck et al. 2 

005 





NGC 4254 

1 

0.5 

- 

Chyzv 

200^ 






NGC 4414 

1 

0.5 

0.5 

iioKla e 

t al. 

l002l 




NGC 6946 

1? 

- 

1? 





.ohde et al. 

1999} 

IC 342 

1 

- 

- 

iiolcolott et a 

.ITT992 


Beclt 

2013 


LMC 

1 

- 

- 

Gaensler et al. |(2U05 

i 

Mao et al.) 

2012 




Fig. 11 Total radio emission (colour) and B-vectors (corrected for Faraday rotation in the foreground 
of the Milky Way) in the And romeda galaxy (M 31), observed at 6cm wavelength with the Effelsberg 
telescope at 3' resolution (from|Giefiubel|j2012)). (Copyright: MPIfR Bonn) 


only due to a regular field in the disk, its line-of sight component should generate 


a conspicuous large-scale pattern in RM, but this is not observed (Fletcher et al. 


20111. This means that a large amount of the ordered field is anisotropic turbulent and 


probably generated by compression and shear of the non-axisymmetric gas flows in 
the density-wave potential. The anisotropic turbulent field is strongest at the positions 
of the prominent dust lanes on the inner edge of the inner gas spiral arms, due to 
compression of turbulent fields in the density-wave shock. A regular held (composed 
of modes m = 0 and m = 2) also exists in the disk of M 51, but is much weaker than 
the anisotropic turbulent held ( Fletcher et al.|[20TT] l. Notably, RMs between 18 cm 
and 22 cm wavelengths, tracing only the polarized emission from regions nearest to 
the observer, show a large-scale pattern that indicates a regular held in the halo of 


M51 dHeald et al.|2009l [Fletcher et al.|201 lt|Mao et al.|2015] l. 
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Fig. 12 Faraday rotation measures between 6 cm and 11 cm wavelengths (colour, from blue=-175 rad m“^ 
to red=+25 rad m“^), total radio emission at 6 cm (contours) and B-vectors (corrected for Fara day rotation) 
in the Andromeda galaxy (M 31), observed with the Effelsberg telescope at 5' resolution (from |Berkhuijsen] 
|et al.lpOO^ ). (Copyright: MPlfR Bonn) 


The large-scale regular field in the nearby galaxy IC 342 is similarly weak as in 
M 51. It is only visible in the pattern of RMs obtained from the low-resolution Effels¬ 
berg images (Beck 2015| l, while RMs derived from the higher-resolution VLA images 
(Fig.[T4|) are about 10 times larger. The field direction changes along the northern spi¬ 
ral arm (Sect. [4T0l i. 

The central regions of M 31 (Fig. 11 1 and IC 342 (Fig. [T7]l host r egular spiral 
fields that are disconnected from the disk fields ( |GieBubel & Beck|2014HBeck|2015| l. 
As the direction of the radial field component points outwards, opposite to that of the 
disk field, two separate dynamos seem to operate in these galaxies. 

In the disks of many other galaxies no clear pattern of Faraday rotation was found. 
Either several high-order modes are superimposed and cannot be distinguished with 
the limited sensitivity and resolution of present-day telescopes, or the timescale for 


the generation of large-scale modes is longer than the galaxy’s lifetime (Arshakian 


et al.||2009]l. Also field injection by strong star-formation activity may perturb the 


generation of a large-scale regular field ( Moss et al.|2012| l. 

While the azimuthal symmetries of regular magnetic fields have been well ob¬ 
served in many galaxies (Table|^, the vertical symmetry (even or odd) is much harder 
to determine. The field of odd-symmetry modes reverses its sign above and below the 
galactic plane. The symmetry type is best visible in strongly inclined galaxies, via 
the RM signs above and below the galaxy plane. In mildly inclined galaxies the RMs 
of diffuse polarized emission from even and odd-symmetry fields differ by a factor 
of about two, which is hard to measure, while the RMs of background sources are 
close to zero for even symmetry (because the RM contributions from the two halves 
cancel) but large for odd symmetry. Background RMs in the area of the Large Mag¬ 
ellanic Cloud (LMC) indicate an even-symmetry field ( |Mao et al.|2012| l. Indications 
for even-symmetry patterns were found in NGC 891 ( |Krause|2009 |l and NGC 5775 
( |Soida et al.|201 1) . 

Spectro-polarimetric data of spiral galaxies are still rare and the application of 
Faraday Synthesis (Sect. |3.6| l has just started (e.g. IGieBiibel et al.|[20T3j |Mao et aT] 
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|2015| ). The Faraday spectrum reflects the distribution of regular magnetic fields, ther¬ 
mal gas and CREs and in principle allows tomography of the ISM in the disk and 
the halo. However, Faraday spectra do not have a simple correspondence to physical 
properties and cannot provide a straight-forward interpretation, as demonstrated by 
the models of|Ideguchi et al. (2014|l. 


4.8 Testing dynamo models 


Observations provide several qualitative evidences for the action of dynamos in spiral 
galaxies: 

(1) The tight radio-lR correlation (Sect. |4.4[ ) and the similarity of kinetic and mag¬ 
netic energy densities (Fig. indicate held amplification by turbulent gas motions 
generated by star-forming processes (|Tabatabaei et al.||2013b[ [Schleicher & Beck| 
[2013] ). 

(2) The magnetic spiral patterns observed in all spiral galaxies so far indicate a gen¬ 
eral decoupling between magnetic fields and the (almost circular) gas flow due to 
magnetic diffusivity. The magnetic pitch angle deviates systematically from that of 
the spiral structures (e.g. in M 83, Fig.[^, as predicted by the a — Q dynamo. Mag¬ 
netic spiral patterns also exist in flocculent galaxies (Sect. 4.14 1 . At present, no other 


model can explain the magnetic spiral patterns in the various types of galaxies. 

(3) Large-scale regular fields are observed in all spiral galaxies so far (Table |^, an 
important result giving a strong hint to the action of a —Q dynamos. No alternative 
model exists so far. 

(4) By comparing the signs of the RM distribution and the velocity held on both sides 
of a galaxy’s major axis, the inward and outward directions of the radial component 
of the axisymmetric spiral held can be distinguished ( [Krause & Be^ 1998| l. Dynamo 
models predict that both signs should have the same probability, which is confirmed 
by observations. The axisymmetric fields of M 31, IC 342, NGC 253, and the ax¬ 
isymmetric held component in NGC 6946 point inwards, while those of NGC 891, 
NGC 4254 and NGC 5775 ( Krause||2009 1, and the axisymmetric component of the 
disk held in M 51 ( Fletcher et al.|201 l| l point outwards. 

A major effort of quantitative comparisons between observable quantities and 
predictions was performed by [Van Eck et al. (2015 i. Identifying the precise form of 
the connection between magnetic and other galactic properties is not straightforward 
because the parameters can combine to produce non-trivial scalings. For example, 
different physical mechanisms have been proposed to provide saturation of the held 
growth in dynamo models, like the balance between magnetic and kinetic (turbu¬ 
lent) forces, the balance between Coriolis and Lorentz forces, or outflows of hot gas 
carrying small-scale magnetic fields out of the disk in order to preserve the balance 
of magnetic helicity. These non-linear saturation mechanisms depend differently on 
galactic parameters, such as star-formation rate, gas density, rotational velocity, and 
rotational shear. 

From their sample of 20 well-observed galaxies. Van Eck et al. ( [2015 1 found a 
statistically significant relation of the total magnetic held strength with the surface 
density of molecular gas and surface density of the star-formation rate, confirming 
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NGC6946 6cm Polarized Intensity B-Vectors (VLA+Effelsberg) -t- H-Alpha 
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Fig. 13 Polarized radio emission (contours) and B-vectors of NCC 6946, combined from observ ations 
at 6 cm wavelength with the VLA and Effelsberg telescopes and smoothed to 15" resolution (from [Beck] 
j2007)), overlaid onto an Ha image from Anne Ferguson. (Copyright: MPIfR Bonn) 


earlier results (Sect. 4.4 1 . The only other significant relation was found between the 
magnetic pitch angle and the strength of the axisymmetric spiral (m = 0) component 
of the regular field; a more tightly wound held has a stronger ASS component. How¬ 
ever, IC 342 does not ht into this relation; has a tightly wound spiral held (Fig. 10 1 , 
but only a weak axisymmetric component ( |Bec k|2015| l. 

In conclusion, dynamos certainly operate in all spiral galaxies, but other processes 
are amplifying and shaping the held as well. MHD density waves are compressing 
and aligning the held along the gaseous spiral arms. Non-axisymmetric gas hows 
around spiral arms and bars are shearing held lines. Parker instabilities form helically 
twisted held loops that are winding around spiral arms. To measure the importance 
of these competing or cooperating effects, the fundamental scaling relations between 
the properties of magnetic helds and the other galactic parameters need to be deter¬ 
mined, based on radio observations of a large number of galaxies with high angular 
resolution. 


4.9 Magnetic arms 

Surprisingly, in most spiral galaxies observed so far the highest polarized intensi¬ 
ties (i.e. the strongest ordered helds) are detected between the optical arms, hlling a 
large fraction of the interarm space, sometimes concentrated in magnetic arms, as in 
NGC 6946 (Fig. [T3] l. The southern-sky spiral galaxy NGC 2997 hosts compressed 
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magnetic fields at the inner edges of material arms, as well as one well-developed 
magnetic arm ( |Han et ah |1999| l. The strong density waves in NGC 2997 lead to 
high degrees of polarization at 6 cm wavelength of typically 25% at the inner edge 
of the northern arm, similar to those in M 51 ( Fletcher et al.|[20TT] l. The degree of 
polarization of the magnetic arm of 40% and its length of at least lOkpc are similar 
to the values in NGC 6946 and are much higher than in IC 342 ( Beck|[20T5] l. Like 
IC 342, M 51 has a short, rudimentary magnetic arm with a low degree of polar¬ 
ization ( Fletcher et al.|20lT I. NGC 1566 shows signatures of a magnetic arm in the 
south-east ( Ehle et al.||1996 1 , but the angular resolution of these radio observations 
was too coarse to clearly detect magnetic arms. In large barred galaxies, ordered 
magnetic helds also fill most of the interarm space, e.g. in NGC 2442 ([Harnett et al. 


|2004| ), NGC 1097, and NGC 1365 ( [Beck et al.|2005| ), but do not form well-defined 
magnetic arms. 

Several mechanisms have been proposed to explain the high degree of held order 
in interarm regions: 

(1) Magnetic held ropes as a result of a magnetic buoyancy instability in a turbulent 
high-j3 plasma (Kleeorin & Rogachevskii 19901; however, the ISM of galaxies is a 
low-)3 plasma, meaning that the thermal pressure is lower than the magnetic pressure 
(Sect.l^; 

(2) slow MHD density waves in the rigidly rotating inner region of a galaxy (|Lou 


et al.|199'9||Poedts & Rogava|2002|l; however, slow MHD density waves in 3D can 


be subject to the Parker and shearing instabilities ( jFoglizzo & Tagger|1995j l; 

(3) more efficient action of the mean-held dynamo between the optical arms due to 
lower turbulent velocity in interarm regions ( |Moss|19W Shukurov|1998| ); however, 
the observed turbulent velocity is not lower in the interarm regions (e.g. in IC 342, 
jCrosthwaite et al.| ( |200I| l); 

(4) introduction of a relaxation time of the magnetic response in the dynamo equation. 


leading to a phase shift between the material and magnetic spiral arms (Chamandy 


jet al.||20I3a|b| l; however, the resulting magnetic arms are restricted to a relatively 
small region around the corotation radius and have a much smaller pitch angle than 
the gaseous arms; 

(5) drift of magnetic helds with respect to the gaseous arms in a non-axisymmetric 
gas How caused by a spiral perturbation (Otmianowska-Mazur et al. 2002| l or by a bar 
( jKulpa-Dybel et al.|201 Ij l; 

(6) weakening of the mean-held dynamo in the material arms by continuous injection 
and amplification of turbulent helds by supernova shock fronts ([Moss et al.||2013 
[20T5l ); 

(7) weakening of the mean-held dynamo in the material arms by star-formation- 
driven outhows ( jChamandy et al.|2015| l; 

(8) magnetic arms as a transient phenomenon during the evolution of galactic mag¬ 
netic helds, possibly related to the short lifetimes of spiral patterns seen in numerical 
simulations (e.g. Dobbs & Baba|20r4 l. 

Models (5) - (8) are promising. Still, these models are simplihed and consider 
either gravitational perturbations or dynamo action. Self-consistent MHD models of 
galaxies including the gravitational potential with spiral perturbations and mean-held 
dynamo action are still missing. 
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The origin of the rudimentary magnetic arm in IC 342 (Fig. deserves a detailed 
discussion. The absence of long magnetic arms in IC 342 is surprising in view of the 
apparent similarities between IC 342 and NGC 6946, like the rotation curve and the 
star-formation surface density Zsfr (Calzetti et al. ]20T0l l. The evolution of a spiral 
magnetic pattern in IC 342 may be hampered by the possible tidal interaction with the 
Local Group ( |Buta & McCall! 1 999| l and/or by the complex spiral pattern in the gas 
of IC 342, described by a superposition of a two- and four-armed spiral pattern with 
different pattern speeds ( |Meidt et al.||2009] l, which may lead to an unstable spiral 
pattern. The mean-field dynamo needs a least a few rotation periods to build up a 
regular field (e.g. Moss et al.|2012| l. A short lifetime of a stable pattern may not allow 
the formation of magnetic arms. 

The spiral pattern in the gas of M 51 is also distorted and short-lived, so that 
only a weak large-scale regular field and a rudimentary magnetic atm could develop. 
On the other hand, NGC 6946 has a two-armed spiral pattern with a well-defined 
pattern speed (Fathi et al.| |2007| l, the spiral field extends smoothly into the central 
region ( |Beck|[^07| l and the large-scale regular field is strong ( Ehle & Beck||1993] l. 
In summary, the existence of magnetic arms may indicate a stable spiral arm pattern 
over several galactic rotation periods, which may not be given in IC 342. 

Most models of spiral gravitational perturbations in stars and gas have neglected 
the effect of magnetic fields so far. On the other hand, the only dynamo model for 
the amplification and ordering of magnetic fields including the spiral perturbations of 
gas density and gas velocity is the kinematical model by|Otmianowska-Mazur et ak 


( 2002| l, in which the back reaction of the field onto the gas flow was not included. The 
MHD model of Pakmor & Springel| ( |201^ includes self-gravity and spiral perturba¬ 
tions, but no mean-held dynamo action. There is urgent need for a synergy between 
these approaches to achieve a comprehensive description of the evolving magnetized 
ISM in galaxies. 


4.10 Helically twisted Parker loops 


IC 342 is one of the best studied nearby spiral galaxies. Its small distance and bright 
radio disk allows detailed studies of the magnetic held that are still impossible in 
most other galaxies. 


Faraday rotation measures (Fig. 14 1 are found to vary periodically along the most 
prominent northern radio spiral arm. The average distance between the extrema cor¬ 
responds to about 2.1 kpc. This feature can be interpreted as a magnetic loop (Parker 
instability) extending out of the galaxy’s disk and out of the sky plane, giving rise to 
a periodic pattern in RM. A regular held bending out of the plane should lead to min¬ 
ima in polarized intensity (tracing the held component in the sky plane) at locations 
where \RM\ (tracing the held component along the line of sight) is at its maximum. 
Minima of polarized intensity occur roughly at the locations of the extrema in RM, 
with a similar distance between the maxima or minima of about 2.2 kpc. 

This is the second indication of a Parker loop in the magnetic held of a nearby 
galaxy, after M 31 ( |Beck et ar]|1989| l. The numerical models by |Kim et al.| ( [2002| ) 
(who assume (5 = 1) predict a wavelength of the most unstable symmetric mode 
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Fig. 14 Faraday rotation measures (RM) between 3.5 cm and 6 cm wavelengths (colour) and total emission 
at 6 cm (contours) in the central and northern regions of IC 342, observed with the VLA at 25 " resolution. 
RM is shown only at points where the polarized intensities at both wavelengths exceeds 3 times the rms 
noise (from |Beck|j2015t ). 


between AnH and 17//, where H is the scale height of the gas. The peak-to-peak 
wavelength of about 4 kpc measured in IC 342 corresponds to // « 230 — 320 pc. 


which is larger than typical scale heights of Hi gas disks of spiral galaxies (Bagetakos 


et al.||20lT I. The reason for this discrepancy may be that the assumption of j3 = 1 i 


IS 


incorrect, because observations indicate j3 < 1 (Sect. 4.2 1 . 

The magnetic field in the northern arm of IC 342 diverts not only in the vertical di¬ 
rection, but also in the disk plane. The ridge line of the polarization spiral arm as well 
as the magnetic pitch angle (Fig. oscillate around the northern arm with similar 
periodicity. This gives indication for a large-scale, helically twisted flux tube, as pre¬ 
dicted by models of the Parker loops under the Coriolis force (jShibata & Matsumoto 
[TWT]|Hanasz et al.|2002| l. 

In the galaxy NGC 7479 (Fig.[24)l, two jets serve as bright polarized background 


sources. High-resolution observations revealed several reversals in RM on 1 — 2 kpc 
scale, originating in the foreground disk of the galaxy ( |Laine & Beck|2008] ), which 
may represent another case of a helically twisted held loop. 


4.11 Large-scale held reversals 

A large-scale held reversal exists between the central region and the disk in M 31 
( [Giefitibel & Beck|2014| l and in IC 342 ( |Beck|[20T5] l. A large-scale held reversal at 
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NGC1097 6.2cm VLA Total Intensity + B-Vectors (PI) HPBW=10'' 
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Fig. 15 Total radio emission (contours) and B-vectors of the barred galaxy NGC 1097, observed at 6 cm 
wavelength with the VLA and smoothed to 10^^ resolution (from [Beck et al.|j2005) ). The background 
optical image is from Halton Arp. 


about constant azimuthal angle across the whole galaxy disk was found in NGC 4414 
( Soida et al.||2QQ2| >. To fulfil the divergence-free condition, this reversal cannot oc¬ 
cur in the same volume, but e.g. between disk and halo. Indication for a large-scale 
field reversal between disk and halo was also found in M 51 (jBerkhuijsen et aL|1997| 
[Fletcher et al.|2Ql l||M^et al.|2Q15| >. New observations and applying Faraday Syn- 
thesis (Sect. |3.6) are needed. 

Large-scale field reversals between spiral arms, like that observed in the Milky 
Way (Fig. [26| , have not been detected in external galaxies, although high-resolution 
images of Faraday rotation are available for many spiral galaxies. Possible reasons 
are discussed in Sect.|5] 


4.12 Barred galaxies 


Galaxies with massive bars have non-axisymmetic gas flows that interact with mag¬ 
netic fields. The magnetic field lines in NGC 1097, one of the closest and brightest 
barred galaxies (Fig. 15 i, follow a pattern that resembles the flow of cold gas in a bar 
potential. As the gas rotates faster than the bar pattern of a galaxy, a shock occurs in 
the cold gas, which has a small sound speed, while the flow of warm, diffuse gas is 
only slightly compressed but sheared. The ordered field is also hardly compressed. 
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NGC10d7-Center 3.5cm VLA Total Int. + B-Vectors (PI) HPBW=3" 



Fig. 16 Total radio emission (contours) and B- 
vectors in the circumnuclear ring of the barred galaxy 
NGC 1097, observed at 3.5 cm wavelength with the 
VLA at y resolution. The background optical image is 
from the HST (from |^ck et al. 1^2005) ). 
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Fig. 17 Polarized intensity and B-vectors in 
the central region of IC 342, observed at 3.5 cm 
wavelength with the VLA at 12" resolution, 
overlaid on a greyscale presentation of the 
CO(l-O) emission at 2.6 mm from the BIMA 
survey (Heifer et al.|2003) . 


probably coupled to the diffuse gas and strong enough to affect its flow (Beck et al. 


20051. The ordered held is also strong in the upstream region (south of the center 


in Fig. 15 I, oriented almost perpendicularly to the bar and aligned with the dust hla- 
ments seen on the optical image in the background. The polarization pattern in barred 
galaxies can be used as a tracer of shearing gas Hows in the sky plane and comple¬ 
ments spectroscopic measurements of radial velocities. 

The central regions of barred galaxies are often sites of ongoing intense star for¬ 
mation and strong magnetic Helds that can affect gas flows. NGC 1097 hosts a bright 
ring with about 1.5 kpc diameter and an active nucleus in its center (Fig. [T6| ). The 
ordered Held in the ring has a spiral pattern and extends to the nucleus. The orienta¬ 
tion of the innermost spiral Held agrees with that of the spiral dust fllaments visible 
on optical images. Magnetic stress in the circumnuclear ring due to the strong total 
magnetic Held (about 50 /rG) can drive gas inflow ( Balbus & Hawley|1998| l at a rate 
of several solar masses per year, which is sufficient to fuel the activity of the nu¬ 
cleus ( |Beck et al.|200 5[). MFID modeling confirmed that magnetic Helds can strongly 
enhance the mass gas inflow rate (Kim & Stone|2012 1. 


The central region of IC 342 hosts a bar of dust and cold gas with a wealth of 
molecular lines. In terms of size, dynamical mass, molecular mass, and star-formation 
rate, the nucleus of IC 342 is a potential twin of the Galactic centre ( Meier|[2()14| l. 
While the central bar in total radio intensity coincides with the central bar in the 
CO line emission from molecular gas, the polarized emission reveals a double-lobe 
structure that is displaced from the CO bar (Fig.[T7]i. 
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Fig. 18 Total radio emission and B-vectors 
of the edge-on galaxy NGC 891, a galaxy 
similar to the Milky Way, observed at 3.6 cm 
wavelength with the Effelsber g telesc ope at 
84'' resolution (from [Krause ] <2009) ). The 
background optical image is from the CFHT. 
(Copyright: MPIfR Bonn) 


Fig. 19 Total radio intensity and B-vectors of the edge-on 
galaxy NGC 5775, combined from observations at 3.6 cm 
wavelength with the VLA and Effelsberg telescopes at 16” 
resolution (from|Soida^etalJJ20UJ)^^ The background Ha 
image is from |TUllmann et al.H2000) . (Copyright: MPIfR 
Bonn) 


4.13 Magnetic halos 


Hot gas, magnetic fields and cosmic rays in galaxy disks drive outflows that extend 
far beyond the galaxy disks seen at optical wavelengths and form thick disks or halos 
at radio wavelengths. Radio polarization observations of several edge-on galax¬ 
ies reveal vertical field components in the halo forming an X-shaped pattern, like 


in NGC 253 ( jHeesen et ai:i|2009b| l, NGC 891 (Fig. [g, NGC 4631 (Fig. [20|, and 
NGC 5775 (Fig. |19|l, which may be related to outflows or to dynamo action ( Moss| 
|et al.|2010| ) driven by such outflows. 

The vertical profile of the total radio continuum emission gives information about 
the outflow speed. The profiles of five edge-on spiral galaxies observed with high res¬ 
olution can be described by two exponential scale heights H^yn, 300 ±50 pc for the 
thin disk and 1.8±0.2kpc for the halo ( Krause|2014| l. Neglecting the thermal con¬ 
tributions and assuming energy equipartition between total magnetic fields and total 
cosmic rays (Sect. |3.4| l, the scale heights of the total magnetic field are typically 
Hb ± (3 ± a) i/syn ~ 1.2 kpc for the thin disk and about 7 kpc for the halo, but proba¬ 
bly larger due to energy losses of CREs. Because the average field strengths and hence 
the synchrotron lifetimes of CREs are different in different galaxies, the roughly con¬ 
stant halo scale heights indicate that outflow speeds increasing with the average field 


^ As a definition of how to distinguish “thick disk” from “halo” is missing, “halo” is used in the follow¬ 


ing. 
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Fig. 20 Total radio emission (contours) and B-vectors of the edge-on g alaxy NGC 4631 , obser ved at 
3.6 cm wavelength with the Effelsberg telescope at 85^^ resolution (from [Mora & Krause| (loT^ ). The 
background optical image is from the Digital Sky Survey. (Copyright: MPIfR Bonn) 
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strength (and hence with the star-formation rate) may balance the smaller CRE life¬ 
times ( |Krause|2009[ ). 

Stronger magnetic fields in the central regions cause higher synchrotron loss of 


CREs, leading to the “dumbbell” shape of many radio halos, e.g. in NGC 253 (Heesen 


et al.|2009all and NGC 4565 (Krause|2009|l. Erom the radio scale heights at several 


frequencies and the corresponding lifetimes of CREs (depending on synchrotron, in¬ 
verse Compton and adiabatic losses) a wind speed of about 300 km/s was measured 
for the electrons in the halo of NGC 253 ( |Heesen et al.|2009a| . 

Galaxies without a detectable radio halo are rare. UGC 10288 shows discrete 
high-latitude radio features, but it does not have a global radio continuum halo El 
win et al. 20131. No radio halo was found in M 31 (|Berkhuijsen et ah 20131 and 


in NGC 7462 (Heesen et al., MNRAS, in press). These three galaxies have a star- 
formation rate that is probably too low to drive outflows. 

On the other hand, exceptionally large and almost spherical radio halos are ob¬ 
served around interacting galaxies, namely NGC 4631 (Eig. [20| and the starburst 
galaxy M 82 (Adeba hr et al.|2013][Reuter et al.|1994 |). A few magnetic spurs could be 
resolved in NGC 4631, connected to star-forming regions ( |Golla & Hummel|1994| l. 
These observations support the ideas of a strong galactic outflow {galactic wind), 
driven by regions of star formation in the disk, possibly supported by the interaction 
with neighbouring galaxies. 

Above the central starburst region of NGC 253, an outwards-directed helical field 
of about 20 flG strength, extending to at least 1 kpc height, could be identified in the 
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gas outflow cone with help of high-resolution RM mapping ( |Heesen et al.||20lTal l. 
This held may help to confine the outflow. 

In summary, the occurrence of radio halos is governed by the competition be¬ 
tween vertical transport of CREs by diffusion and/or convection (galactic wind) and 
energy losses of CREs, dominated by synchrotron loss. CRE transport can be charac¬ 
terized by the escape time fesc needed to reach the scale height H of the disk, beyond 
which the synchrotron emission significantly weakens. In case of diffusive transport, 
fesc where D is the diffusion coefficient, while in case of convective trans¬ 

port fesc ^f/vwindj where Vwind is the average wind speed. 

Eor strong magnetic fields or slow winds or slow diffusion (fsyn fesc), the CREs 
cannot leave the disk and no radio halo develops, as in the case of M 31 ( |Berkhuijsen 
et al.|201^ . Eor tsyn ^ fesc, the scale height // of the synchrotron halo increases with 
distance from the center according to the increase in synchrotron lifetime (“dumb¬ 
bell” halo) and with radio wavelength (// oc for diffusive or H oc for con¬ 
vective propagation). Eor weak fields or fast winds (tsyn > fesc), a large radio halo of 
elliptical or spherical shape is formed with a wavelength-independent scale height 
that is determined by the scale height of the total magnetic field and adiabatic losses 
of the expanding CRE flow.|^ 


4.14 Rudimentary spirals: flocculent and irregular galaxies 


Elocculent galaxies have disks but no gaseous spiral arms. Nevertheless, magnetic 
spiral patterns have been observed in all flocculent galaxies so far, indicating that the 
a — dynamo works independently of spiral density waves. Ordered magnetic fields 
with strengths similar to those in grand-design spiral galaxies have been detected in 
the flocculent galaxies M 33 ( |Tabatabaei et al. 2008|l, NGC 3521 and NGC 5055 
( Knapik et al.|2000l l, and in NGC 4414 ( Soida et al.|2002| l. The mean degree of po¬ 
larization (corrected for the differences in spatial resolution) is also similar in grand- 
design and flocculent galaxies. 

Radio continuum images of irregular, slowly rotating galaxies may reveal strong 


total magnetic fields, e.g. in the Magellanic-type galaxy NGC 4449 (Eig. 21 1 , with 
a partly ordered field of about 7 /rG strength, a spiral pattern in the northeast and a 
radial pattern in the central region ( |Chyzy et al.| 2000 1. Earaday rotation data shows 
that this ordered field is partly regular and the a — dynamo is operating (jSiejkowski 
let al.|2014D . 

The Large and the Small Magellanic Clouds are the closest irregular galaxies at 
distances of 50 kpc and 60kpc, respectively, interacting with the Milky Way and each 
other. Radio observations found the two galaxies to be weakly polarized (|Haynes 
etal.|1991|[Maoetal.|2012| l, with the exception of two highly polarized filaments in 
the southeast of the EMC. A study based on RM grids suggests that both galaxies host 
large-scale coherent fields, indicating that a large-scale dynamo also works under the 
less favourable conditions of slow ordered rotation ( |Mao et al.|2008|[2012| l. Due to 
their low potential well, low-mass galaxies are prone to outflows and galactic winds 
driven by star formation. 


® Detailed models of radio halos are presented in Lerche & Schlickeiser 1 1982a b 
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NGC4449 3.6cm VLA+Effelsberg Tot.Int.+B-Vectors-t-H-alpha HPBW=12” 



Fig. 21 Total radio emission (contours) and B-vectors of the dwarf irregular galax y NGC 4449, observed 
at 3.6cm wavelength with the VLA at 12'^ resolution (from j^yzy et aL]j200Q^ ). The background Ho; 
image is from Dominik Bomans (Bochum University). 


In dwarf irregular galaxies, the strengths of the total (turbulent) field are gener¬ 
ally smaller than in spiral galaxies ( Chyzy et al.|2Qll|>, except for starb urst dwarfs, 
e.g. NGC 1569 with 10 — 15 fJ.G total field strength ( fepley et al.|2QlQ] >, where star 
formation activity is sufficiently high for the operation of the small-scale dynamo. 
Consequently, dwarf galaxies form the low-luminosity tail of the radio-lR correla¬ 
tion for large galaxies ( Chyzy et al.|201 1 1. 

Spots of faint polarized emission were detected in dwarf irregular galaxies, but no 
large-scale regular fields so far (Heesen et al. 201 lb| l, either due to a lack of telescope 
sensitivity or a lack of dynamo action. In classical a — Q dynamo models, the dy¬ 
namo number is too low for dynamo action because rotation is almost chaotic, while 
models including support of dynamo action by outflows predict a—£2 dynamos even 
in dwarf galaxies (Rodrigues et al.|2015|l. 


4.15 Beyond spirals: interacting galaxies 

Gravitational interaction between galaxies leads to asymmetric gas flows, compres¬ 
sion, shear, enhanced turbulence and outflows that can be traced by observing mag¬ 
netic fields in radio polarization. Magnetic fields can become aligned along the com¬ 
pression front or perpendicular to the velocity gradients. Such gas flows make turbu- 
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NGC4535 Effelsberg 3.6cm Pol. Intensity -f B-Vectors (%) HPBW=90" 



Fig. 22 Polarized radio emission (contours) and B-vectors of the Vir go galaxy NGC 4535, ob served at 
3.6 cm wavelength with the Effelsberg telescope at 90^^ resolution (from [^fezgowiec et al.|j2QT^ and priv. 
comm.). The background optical image is from the Digital Sky Survey. 


lent fields highly anisotropic. Large-scale dynamos can be enhanced by ram pressure 
( |Moss et al.|2014| ). 


The classical interacting galaxy pair is NGC 4038/39, the “Antennae” (Chyzy & 


|Beck|20()4| ). It shows bright, extended radio emission filling the whole system. In the 
interaction region between the galaxies, where star formation did not yet start, and at 
the northeastern edge of the system, the magnetic field is partly ordered, probably the 
result of compression and shearing motions along the tidal tail. Particularly strong, 
almost unpolarized emission comes from a region of violent star formation, hidden 
in dust. The average total magnetic field is stronger than in normal spirals, but the 
mean degree of polarization is unusually low, implying that the fields are tangled. The 
Antennae galaxies have been modeled by Kotarba et al. ( 2010| l, who could reproduce 
the main features observed in the radio image. 

The total magnetic field in a sample of 16 systems in different interaction stages 
increases with advancing interaction, which indicates enhanced production of random 
magnetic fields ( [Drzazga et al.|20lT] l. The strength of the ordered magnetic field are 
sensitive tools for revealing global galactic distortions. The pattern of the ordered 
field traces the orientation of gas flow in tidally stretched spiral arms and in tidal 
tails. Such outflows may have contributed to the magnetization of the intergalactic 
medium in the early cosmological epoch. 

Interaction with a dense intergalactic medium also imprints unique signatures 
onto magnetic fields and thus onto the radio emission. The Virgo cluster is a loca¬ 
tion of especially strong interaction effects (Fig. [2^ , and almost all cluster galax¬ 
ies observed so far show asymmetries of their polarized emission because the outer 
magnetic fields were compressed (|Vollmer et al.|2007l|2013[|Wezgowiec et al.|2007| 
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Right AtMnalon (J2000.0) 


Fig. 23 Total radio intensity (contours) and B-vectors 
of the jets in NGC 4258, observed at 3.5 cm wavelength 
with the VLA at 14^' resolution. The background Ha 
image is from the Hoher List Observatory of the Uni¬ 
versity of Bonn (from |Krause & Lohr|j2004^ ). 


Fig. 24 Total radio intensity (contours) and B- 
vectors of the jets in NGC 7479, observed at 
3.5 cm wavelength with the VLA at 8^' reso¬ 
lution. The background shows a Spitzer/IRAC 
image at 3.6/im (NASA/JPL-Caltech/Seppo 
Laine) (from [Laine & Be^j2Q08t . 


|2012[ ). Ordered fields are an excellent tracer of past interactions between galaxies or 
with the intergalactic medium. 


4.16 Spiral galaxies with jets 


Nuclear jets are observed in only a few spiral galaxies. These jets are weak and small 
compared to those of radio galaxies and quasars. Detection is further hampered by the 
fact that they emerge at some angle with respect to the disk, so that little interaction 
with the ISM occurs. If the nuclear disk is oriented at a large angle to the disk, the 
jet hits a large amount of ISM matter, CREs are accelerated in shocks, and the jet 
becomes radio-bright. However, not all jets are consistent with this morphology. 

NGC 4258 is one of the rare cases where large radio jets of at least 15 kpc length 


are observed (Fig. 23 i. A nuclear disk is observed in water maser emission, has an 


inner radius of 0.13 pc and is seen almost edge-on ( Greenhill et al.|[T995| l. As the 
jets emerge perpendicular to the nuclear disk, they have to pass the galactic disk at 
a rather small angle. The magnetic field orientation is mainly along the jet direction. 
The equipartition field strength is about 300/iG at the resolution of about lOOpc 
( Krause & L6hr|2004| l, which is a lower limit due to energy losses of the CREs and 
the limited resolution. 

The barred galaxy NGC 7479 also shows remarkable jet-like radio continuum fea¬ 
tures: bright, narrow, 12 kpc long in projection, and containing an aligned magnetic 
field (Fig. [24|). The lack of any optical or near-infrared emission associated with the 
jets suggests that at least the outer parts of the jets are extraplanar features, although 
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close to the disk plane. The equipartition strength is 35 — 40 flG for the total magnetic 
field and about 10/rG for the ordered magnetic field in the jets. According to Fara¬ 
day rotation measurements, the large-scale regular magnetic field along the bar points 
towards the nucleus on both sides. Multiple reversals on scales of 1 — 2kpc are de¬ 
tected, probably occurring in the galaxy disk in front of the eastern jet by anisotropic 
fields in the shearing gas flow in the bar potential. 

Highly polarized radio emission from kpc-size jets has also been detected in 
NGC 3079 ( [Cecil et^|2001[ ), with the held orientations perpendicular to the jet’s 
axis, and in the outflow lobes of the Circinus Galaxy ( Elmouttie et al.|[l995| l. Two 
members of the Virgo cluster, NGC 4388 and NGC 4438, have elongated radio fea¬ 
tures emerging from the active Seyfert-type nuclei and extending roughly perpendic¬ 
ular to the planes of the disks ( Hummel & Saikia|1991| l. Detailed images including 
polarization are forthcoming as results of the CHANG-ES survey of edge-on galaxies 
( jlrwin et al.|2012| ). 

Jets in spiral galaxies may be more frequent than the available radio observations 
suggest. Euture low-frequency observations will help to search for weak synchrotron 
emission from interface regions between the jets and the low-density halo gas. 


5 How does the Milky Way fit into the picture of nearhy galaxies? 


Direct measurements of the Voyager 2 spacecraft in the heliosheath indicate that the 
surrounding interstellar magnetic held is 4 — 5 pG strong and oriented at an angle 
of about 30° from the Galactic plane ( Opher et al.|2009l l, probably because the ISM 
held twists close to the heliosphere. Voyager 1 crossed into interstellar space in 2012 


and measured a smooth increase in held strength to 5.62 ± 0.01 pG (Burlaga et al. 


|2013| ). This value is very close to those obtained with other methods (see below) and 
to those in nearby galaxies with a low star-formation rate (Sect. 4.1 1 . 

To deduce information about the strength and structure of the Galactic magnetic 
held, pulsars are ideal objects because their RMs provide held directions at many dis¬ 
tances from the Sun ( |Noutsos|2012 1. Since most pulsars are concentrated along the 
Galactic plane, they sample the held in the disk. Combination of RM and DM data 
of pulsars (Eq. assuming uncorrelated fluctuations, gives an average strength of 
the local regular held along the line of sight of 2.1 ± 0.3 pG and about 4 pG at 3 kpc 
Galacto-centric radius ( Han et al.|2006 l. These are upper limits in the case of corre¬ 
lated fluctuations ( [Beck et al.|2003| l,while lower limits in the case of anticorrelated 
fluctuations or field reversals along the line of sight. 

Erom the dispersion of pulsar RMs, the Galactic magnetic held was found to have 


a significant turbulent component with a mean strength of 5 — 6 pG (Rand & Kulkarni 


|1989[ jOhno & Shibata|[l993[ [Han et ar]|2004| ). Magnetic turbulence occurs over a 
large spectrum of scales, with the largest scale determined from pulsar RMs of Iturb — 
55 pc ( [Rand & Kulkarni| 1989) 1 or — 10 — 100pc ( jOhno & Shibata|1993| l. These 
values are consistent with the size of turbulent cells of ~ 50 pc estimated from beam 


depolarization and Earaday depolarization in external galaxies (Sects. 3.5 and 3.7 1 . 
However, turbulence scales of only a few parsecs in spiral arms were derived from 
RM structure functions of polarized background sources ([Haverkorn et al.|2008jl. 
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Fig. 25 All-sky map of rotation measures in the Milky Way, constructed from the RM data of about 40000 
polarized extragalactic sources from the VLA NVSS sur vey jXaylor et al.|2009^ and other catalogs. Red: 
positive RM, blue: negative RM (from[Oppermann et al.H2012^). 


Modeling the surveys of the total synchrotron and y-ray emission from the Milky 
Way yield field strengths near the Sun of about 5 fiG of the isotropic turbulent field. 


about 2 /iG of the anisotropic turbulent field and about 2 /rG of the regular field (Or- 


lando & Strong|20T3)l, adding up to a total field strength of about 6 /rG. This is in ex¬ 


cellent agreement with the Voyager and pulsar RM data (see above) and the Zeeman 
splitting data of low-density gas clouds ( [Crutcher et al.||2OT^ . In the inner Galaxy 
the total field strength is about 10 /rG. In the synchrotron filaments near the Galactic 
Center the total field strength is about 100/dG ( Crocker et al.|201^ . 

Optical polarization data of about 5500 selected stars in the Milky Way yielded 
the orientation of the large-scale magnetic field near the Sun ( Fosalba et al.||20()^ , 
which is mostly parallel to the Galactic plane and oriented along the local spiral arm. 

The all-sky maps of polarized synchrotron emission at 1.4 GHz from the Milky 
Way from DRAG and Villa Elisa and at 22.8 GHz from WMAP, and the Effelsberg 
RM survey of polarized extragalactic sources, were used to model the regular Galactic 
field ( jSun et al.|2008| Sun & Reich|2010| l. One large-scale reversal is required at 
about 1—2 kpc from the Sun towards the Milky Way’s center, consistent with pulsar 
data. 


RM data from pulsars and extragalactic radio sources (Eig. 251 was used to model 


the Galactic magnetic field ( Nota & Katgert|2010 Van Eck et al.||20lT] Jansson & 


|Farrar|2012D . A large-scale magnetic field reversal appears to be present between the 
Scutum-Crux-Sagittarius arm and the Carina-Orion arm (Fig. 26 1 . As distances to 
most pulsars are uncertain, this result should be taken with some caution. The overall 


structure of the regular field in the disk of the Milky Way is still uncertain (Noutsos 


20091. A larger sample of pulsar RM data and improved distance measurements to 


pulsars are needed. A satisfying explanation for the large-scale reversal in the Milky 
Way is still lacking. So far no similar reversals have been detected in external galaxies 
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Fig. 26 Model of the magnetic field in the Milky Way (vectors), derived from Faraday rotation measures of 
pulsars and extragalactic sources. Generally accepted results are indicated by yellow vectors, while white 
vectors are not fully confirmed. The background image shows the distribution of ionized gas delineating 
the spiral arms (from Jo-Anne Brown, Calgary). 


(Sect. [4TT] i. Possible reasons are: 

(1) The reversal in the Milky Way may be restricted to a thin layer near to the plane 
and therefore hardly visible in the average RM data of external galaxies along the line 
of sight. 

(2) The reversal in the Milky Way may be of limited azimuthal extent and difficult to 
observe in external galaxies with present-day telescopes. 

(3) The reversal in the Milky Way may be part of a disturbed field structure, e.g. 
due to interaction with the Magellanic clouds, or a relic from seed fields ( |Moss et ak] 

[M2l ). 


Even less is known about the halo field of the Milky Way. The vertical full equiv¬ 
alent thickness of the synchrotron emission from the Milky Way is about 3 kpc near 
the Sun ( |Beuermann et al.| 19851 1, scaled to a distance to the Galactic center of 8.5 kpc, 
yielding an exponential scale height of //syn — 1 5 kpc and at least 6 kpc of the total 
field, similar to that in external galaxies (Sect. 4.13| l. 

The signs of RMs of extragalactic sources and of pulsars at Galactic longitudes 
I — 90° — 270° are mostly the same above and below the plane (Fig. [25] l: the mag¬ 
netic field in the local disk is symmetric, while the RM signs towards the inner Galaxy 
{I = 210° — 90°) are opposite above and below the plane. This could be assigned to an 
odd-symmetry (antisymmetric) halo field ( Sun & Reich|2010) l. However, according 
to RM data from extragalactic sources, the local regular Galactic field has no signifi¬ 
cant vertical component towards the northern Galactic pole and only a weak vertical 
component of 5^ ~ 0.3 /rG towards the south (Mao et al. 2010| l. This is neither consis¬ 
tent with an odd-symmetry nor with an even-symmetry halo field as found in several 
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Table 3 General field structures in nearby galaxies 


Galaxy type 

Magnetic field structure 

Regular Held 

Sc with strong 

Spiral at inner arm edge and in 

Strong, 

density wave 

interarm regions, turbulent in arms 

fluctuating 

Sb or Sc with weak 

Spiral in interarm regions, 

Moderate, 

or moderate density wave 

turbulent + ordered in arms 

large-scale 

Barred Sc 

Ordered + turbulent along bar, 
spiral outside bar 

Moderate 

Sa 

Ordered + turbulent 

Detected 

SO 

Not detected 

Not detected 

Interacting spiral 

Ordered, asymmetric 

Weak 

Spiral with nuclear jets 

Ordered along jet 

7 

Flocculent Sc or Sd 

Spiral + turbulent in disk 

Weak 

Irregulai* 

Turbulent in star-forming regions 
- 1 - segments of ordered field 

Weak 

Starburst dwarf 

Turbulent in star-forming regions 

Not detected 

Spheroidal dwarf 

Not detected 

Not detected 

E without active nucleus 

Not detected 

Not detected 


external galaxies (Sect. p^ . The halo field of the Milky Way has a more complicated 
structure than predicted by dynamo models. Modeling the diffuse polarized emission 
and RMs gave evidence for an X-shaped vertical held component pans son & Farrar 


20121, similar to that in external galaxies (Sect. 4.13 i. 


In summary, strength, spiral pattern, and vertical extent of the Galactic magnetic 
held are similar to those in nearby spiral galaxies. Two major differences still need to 
be understood; the large-scale held reversal(s) and the antisymmetric pattern of the 
halo held, both of which are not observed in external galaxies. Either our Milky Way 
is special, or the different observational methods deliver apparently incompatible re¬ 
sults. 


6 Summary and open questions 

Table gives a concise description of the general structures of ordered magnetic 
helds in various galaxy types as observed from radio polarization observations. For 
most galaxy types only a few cases were observed so far. - “Turbulent” stands for 
“isotropic turbulent or tangled”, “ordered” for “anisotropic turbulent” or “regular”. 
“Regular held” refers to additional measurements of Faraday rotation. 


For a list of radio polarization observations of galaxies, see Beck & Wielebinski 


( 2013| l and updated versions on arXiv. Radio images of total intensity and polariza¬ 
tion are compiled on: http://www.mpifr-bonn.mpg.de/atlasmag 


The open questions raised in this review may inspire future work: 


(1) Morphology; 
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- How far do turbulent and regular magnetic fields extend radially in the disk and 
vertically above the disk? 

- How well are magnetic pitch angles aligned with the orientations of gaseous spiral 
arms? 

- How do magnetic pitch angles vary radially in the disk? 

- Do dominating bisymmetric spiral patterns (azimuthal mode m = 1) of the regular 
held exist in galaxies? 

- Do azimuthal modes m>2 exist and how strong are these? 

- Which galaxy types host magnetic arms between gaseous spiral arms? 

- How common are held loops generated by the Parker instability? 

- Do regular helds in galaxies have large-scale reversals? 

- Are the X-shaped helds in galaxy halos related to outhows (anisotropic turbulent 
helds) or do they show large-scale patterns in Faraday rotation (regular helds)? 


(2) Origin and amplihcation; 


- What was the main source of seed helds needed for the amplihcation of galactic 
helds in young galaxies? 

- Are small-scale helds generated by the small-scale dynamo or by tangling of 
regular helds? 

- What is the origin of magnetic arms between gaseous spiral arms? 

- What determines the efficiency and saturation level of the a Q dynamo? 

- Is the mean helicity of the large-scale magnetic helds in galaxy disks non-zero, 
as predicted by a — dynamo action? 

- Are small-scale helds with non-zero mean helicity expelled from the disk to keep 
the a — Q dynamo alive? 


(3) Dynamical importance: 


- Where is energy equipartition between total magnetic held and total cosmic rays 
valid and where does it fail? 

- What is the dominant propagation mechanism of CREs in galaxy disks and halos 
and how does this depend on magnetic held strength and structure? 

- How do magnetic helds interact with spiral density waves? 

- How does the relation between total magnetic held strength and gas density vary 
between galaxies and within galaxies? 

- How do magnetic helds modify the star-formation rate and the star-formation 
efficiency? 

- Is amplihcation of small-scale helds in star-forming regions the reason for the 
radio-IR correlation? 

- Is the magnetic energy density larger than that of gas turbulence in the outer parts 
of galaxies? 

- Can magnetic helds affect the general rotation of gas? 

- Can galactic outhows magnetize the intergalactic medium? 
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Fig. 27 Total synchrotron emission at 20 cm wavelength as a function of redshift z and total magnetic field 
strength B, and the 5 O' detection limits for 10 h and 100 h integration time with the full SKA (from |Murphy| 
^2009) and priv. comm.). 


7 Outlook 


Future high-resolution, high-sensitivity observations at high radio frequencies with 
the Jansky Very Large Array (VLA) and the planned Square Kilometre Array (SKA, 
see below), combined with high-resolution CO observations with the Atacama Large 
Millimeter Array (ALMA), will directly map the detailed field structure and the inter¬ 
action with the molecular gas in external galaxies ( |Beck et al.|2015 1. Images of dust 
polarization of nearby galaxies in the submm range with ALMA will soon become 
possible. Another exciting prospect is extragalactic starlight and infrared polarimetry 
with the European Extremely Large Telescope (E-ELT) ( Strassmeier & Ilyin|2009| l. 

Construction of the first phase of the SKA is planned to start in 2018. The low- 


frequency radio telescopes Low Frequency Array (LOFAR) (van Haarlem et al. 20131 
and the Murchison Widefield Array (MWA) ( |Tingay et al.|2013| are suitable instru- 
ments to search for extended synchrotron radiation in outer galaxy disks and halos 
and the transition region to intergalactic space. Several galaxies have already been 
observed with LOFAR ( Beck et al.|2013[ Mulcahy et al.|2014| ). Low frequencies are 
also suited to search for small Faraday rotation measures from weak fields in the local 
ISM of the Milky Way, in the ISM of external galaxies and in intergalactic space, if 
Faraday depolarization is small. 

Small-scale and large-scale magnetic fields may exist in SO and elliptical galaxies 
without star formation and without active nucleus because turbulence can be gener¬ 
ated by the magneto-rotational instability (MRI) ( Sellwood & Balbus|1999| l or may 
result from star formation activity in the past. The detection may become possible 
with the SKA via deep imaging or RM grids of background sources. 
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The detection of strong radio emission in distant galaxies (which is at least partly 
of synchrotron origin) demonstrates that strong magnetic fields existed already in 
young galaxies with strengths of several 100/iG ( Murphy |2009| l. Total synchrotron 
emission can be detected with the SKA out to even larger redshifts in starburst galax¬ 
ies, depending on luminosity and magnetic field strength (Fig. 27 1 . As the amplifica¬ 
tion of turbulent fields is fast and efficient (Sect.j^ and the star-formation rates were 
high at high redshifts, the total fields can be expected to be very strong already in 
young galaxies. Then the energy range of the electrons emitting at a fixed observa¬ 
tion frequency shifts to low energies (Eq. |^, where ionization and bremsstrahlung 
losses may become dominant, so that there must be a maximum redshift until syn¬ 
chrotron emission is detectable. Its measurement with the SKA will constrain models 
of the evolution of magnetic fields in young galaxies ( [Schleicher & Beck|2013| ). 

Nearby galaxies seen edge-on generally show an ordered, disk-parallel field near 
the disk plane ( Dumke et al.|1995] l. As a result, polarized emission can also be de¬ 
tected from distant, unresolved galaxies if the inclination is larger than about 20° ( jStilj 
jet al.|20()9l l. This opens a new method to search for ordered fields in distant galaxies 
with help of deep surveys. 

If polarized emission from galaxies is too weak to be detected, the method of RM 
grids towards background QSOs can still be applied. Here, the distance limit is given 
by the polarized flux of the background QSO, which can be much higher than that 
of an intervening galaxy that is identified optically by its Mg II absorption system. 
With this method, significant regular fields of several /iG strength were detected in 
distant galaxies ( Bemet et al.|2008 2013| Fames et al.||2014 l. Mean-field solutions 
of the a — Q. dynamo theory predict patches of regular fields observable as RMs in 
evolving regular fields already at z < 3, with increasing coherence scale, until fully 
developed regular fields are formed at z ~ 1 — 1.5 (Arshakian et al.|2009| Rodrigues 


et al.|2015] l. Note that RM values measured today have been reduced due to redshift 
by the dilution factor of (1 -|-z)^^, so that observations at longer wavelengths are 
needed to detect RMs from regular fields at high redshifts. 

Faraday rotation in the direction of polarized lobes of radio galaxies allows to de¬ 


termine the strength and pattern of a regular held in an intervening galaxy (Kronberg 


et al.|199^ Irwin et al.|2013|). RM data indicate that ionized gas may extend to sev¬ 


eral lOkpc distance from galaxies ( jBernet et al.|2013) . A reliable model for the held 
structure in and around galaxies needs RM values from a large number of polarized 
background sources, hence large sensitivity and high survey speed. 

The SKA “Cosmic Magnetism” Key Science Project plans to observe a polariza¬ 
tion survey over the entire accessible sky with the SKA-MID Band 2 (around 20 cm 
wavelength) ( Johnston-Hollitt et al.||2015] l. Within 1 h integration per held this will 
allow the detection of about 10 million discrete extragalactic sources and the mea¬ 
surement of their RMs, about 300 RMs per square degree or a mean spacing between 
sources of about 3'. 

Deep polarization observations with the SKA centered on a sample of nearby 
galaxies will allow studies of the interaction between magnetic fields and the ISM 
gas ( Beck et al.|2015| l. The detailed reconstruction of the 3D magnetic held structure 
will become possible with help of Faraday tomography (Heald et al. |2015j l or a dense 
grid of RM values from background sources ([Stepanov et al.|2d08ll. More than 1000 
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RM values can be expected in the area of M 31. Applying this method to more distant 
galaxies, simple patterns of regular fields can be recognized out to distances of about 
100 Mpc ( Stepanov et al.|200^ . 


Exciting prospects! 
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